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ABSTRACT 


v  — 

A  detailed  investigation  of  the  spatially  anisotropic  com¬ 
ponent  of  the  laterally  averaged  molecular  hydrogen/Ag( 111)  phy- 
sisorption  potential  is  presented.  Experimentally  derived  rota- 
tionally  inelastic  transition  probabilities  for  h£\  d!/,  and  HD, 
taken  as  a  function  of  collision  energy,  are  compared  with  those 
resulting  from  close-coupled  quantum  scattering  calculations. 

These  calculations  utilize  exponential- 3  and  variable  exponent 
parameterizations  of  the  laterally  averaged  isotropic  potential 
which  reproduce  the  experimental  bound  state  resonance  spectra 
for  p-H^  and  o-D^  on  Ag(lll) .  Complementary  information  is 
obtained  by  analyzing  the  magnetic  sub-level  splittings  for  phy- 
sisorbed  J  =  1  n-H^,  using  diffractive  selective  adsorption  re¬ 
sonance  energies  calculated  with  first  order  perturbation  theory. 
Theoretical  predictions  for  HD/Ag(lll)  rotationally  mediated  selec 
tive  adsorption  resonances  are  also  compared  with  previously  re¬ 
ported  experimental  results,  which  show  well  resolved  J-dependent 
energy  shifts  resulting  in  part  from  the  orientational  anisotropy 
of  the  potential.  The  results  obtained  in  this  study  indicate 
that  both  the  attractive  and  repulsive  parts  of  the  anisotropic 
potential  exhibit  only  a  weak  orientation  dependence,  in  agree¬ 
ment  with  recent  theoretical  predictions  for  this  system.  C)Ti  c^/iy 


I.  INTRODUCTION 


Until  recently,  the  phys isorptive  interaction  of  molecules 
with  surfaces  was  far  less  extensively  investigated  than  that  of 
atoms  with  surfaces.  Diffractive  selective  adsorption  experiments 
which  have  been  carried  out  with  molecules,  such  as  those  involving 
H2  and  D2  on  LiF,  NiO,  graphite,  and  stepped  surfaces  of  copper  [1] 
were  typically  analyzed  to  obtain  hound  states  of  the  laterally 
averaged  interaction  potential  with  the  molecule  treated  as  a  struc 
tureless  particle.  Recent  advances  in  molecular  beam  surface  scat¬ 
tering  experimentation  have  now  prompted  a  growing  interest  in  the 
role  of  the  internal  molecular  degrees  of  freedom  in  molecule-sur¬ 
face  interactions.  This  has  been  evidenced  by  a  large  number  of 
experimental  and  theoretical  studies  of  rotationally  and  vibration- 
ally  inelastic  molecule-surface  collision  systems  [2]. 

A  prime  objective  in  these  studies  is  the  determination  of 
the  molecule-surface  interaction  potential.  Because  of  the  mole¬ 
cular  orientation  dependence  this  is  considerably  more  complicated 
for  a  molecule-surface  collision  than  for  atom-surface  collisions. 
For  a  diatomic  molecule  incident  upon  a  perfectly  periodic  surface, 
the  molecule-surface  interaction  potential  may  be  written  as  the 
Fourier  sum 

V(z,  6,  ♦,  R,  r)  *  L  VG(.z,  6,  <*>,  r)  e1?’?  (1.1) 

where  z  is  the  normal  distance  from  the  molecule  center  of  mass  to 
the  surface,  R  ■  (x,  y)  is  the  surface  projection  of  the  center  of 
mass  relative  to  some  origin  on  the  surface,  0  and  <f>  are  the  polar 
and  azimuthal  orientation  angles  of  the  molecular  orientation  in 


the  xyz  frame,  and  r  is  the  molecular  internuclear  distance 
(Figure  1) .  We  shall  be  concerned  with  low  energy  (<  150meV) 
collisions  of  the  molecular  hydrogen  isotopes  for  which  a  rigid 
rotor  assumption  is  valid.  The  laterally  averaged  potential 
Vqq(z,  0,  <j>,  r)  can  then  only  be  a  function  of  z  and  6,  since 
physically  there  can  be  no  azimuthal  orientation  dependence  for 
a  diatomic  interacting  with  a  flat  surface.  Changes  in  rotational 
polarization  or  rotational  state  selected  molecules  can  only  arise 
from  the  presence  of  higher  order  diffractive  terms  (G  +  0)  in 
Equation  (1.1).  (However,  matrix  elements  of  the  laterally  averaged 
potential  Vq0(z,  ®)  will  be  different  for  different  magnetic  quantum 
states  of  the  diatom,  with  the  result  that  an  overall  apparent  change 
in  rotational  polarization  may  be  produced  by  scattering  an  initial 
ensemble  of  | Jm>  molecules  from  a  flat  surface  [3].)  We  restrict 
our  attention  here  to  the  laterally  averaged  potential  VQ0(z,  9): 
this  is  the  component  probed  by  rotationally  inelastic  scattering 
from  smooth  surfaces  and  by  selective  adsorption  resonance  pheno¬ 
mena  on  weakly  corrugated  surfaces. 

Several  electron  energy  loss  spectroscopy  (EELS)  and  molecular 
beam  scattering  experiments  recently  carried  out  for  the  adsorption 
of  molecular  hydrogen  on  smooth  metal  surfaces  indicate  that  the 
physisorbed  molecule  behaves  as  a  nearly  free  3-dimensional  rotor, 
with  only  weak  perturbations  due  to  the  anisotropy  of  the  (laterally 
averaged)  molecule-surface  interaction  [4,  5],  However,  the  low 
energy  resolution  currently  available  in  EELS  experiments  (>  3meV) 
limits  the  capability  of  this  technique  for  investigating  the  exact 
nature  of  the  physisorbed  state  at  the  present  time.  We  have  recently 
reported  a  detailed  comprehensive  study  of  the  diffractive  selective 


3 


adsorption  (hereafter,  DSA)  in  the  molecular  beam  scattering  of 
n-H2»  P‘H2»  n_^2  an<*  °’^2  ^rom  Ag(lll)  which  shows  evidence  of 
small  |J,m>  -dependent  shifts  in  the  bound  states  of  Vqq(z,  6) > 
arising  from  the  weak  anisotropic  component  of  this  potential  [6], 
This  was  made  possible  by  the  high  energy  resolution  of  these 
molecular  beam  experiments  ('v  200yeV)  and  the  development  of  a 
pure  J  =  0  hydrogen  beam.  Such  shifts  have  also  recently  been 
observed  for  on  the  more  strongly  corrugated  Ag(110)  surface 

[7].  Similar  J-dependent  shifts  are  apparent  in  the  recent  measure¬ 
ments  of  rotationally  mediated  selective  adsorption  (henceforth, 
RMSA)  resonances,  for  HD  on  Pt(lll)  [8]  and  for  HD  on  Ag(lll)  [9], 

These  selective  adsorption  studies  are  part  of  a  detailed 
experimental  and  theoretical  study  of  the  molecular  hydrogen/ 

Ag(lll)  physisorption  interaction  potential,  including  the  laterally 
averaged  isotropic*  the  laterally  averaged  anisotropic  and  the 
higher  order  periodic  components.  The  laterally  averaged  isotropic 
component  has  been  treated  in  previous  publications  [6,  10]  while 
the  periodic  potential  will  be  presented  elsewhere  [11].  It  is 
the  purpose  of  this  paper  to  investigate  the  anisotropic  component 
of  VQ0(z,  0)  in  considerable  detail  for  the  molecular  hydrogen 
isotopes/Ag(lll)  system,  employing  a  combination  of  experimental 
measurements  and  theoretical  calculations. 

A  comprehensive  study  of  the  molecular  beam  scattering  of 
hydrogen  from  Ag(lll)  carried  out  in  our  laboratory  has  yielded 
the  following  array  of  information:  rotationally  inelastic  and 
elastic  transition  probabilities  for  H2»  D2,  and  HD,  RMSA  resonances 
for  HD  and  isolated  instances  of  RMSA  resonances  for  H2  and  D2, 
and  DSA  resonances  for  H2  and  D2  in  which  |J,m>  -dependent  shifts 


and  linewidth  differences  (for  are  seen  between  n-F^  and 

an<*  n"^2  anc*  0’^2*  respectively.  Transition  probabilities 
for  simultaneous  diffractive  and  rotationally  inelastic  excita¬ 
tion  of  HD, and  for  diffraction  and  diffraction-rotation  excita¬ 
tion  of  H2  and  D2  have  also  been  measured. 

From  an  experimental  point  of  view  the  low  energy  scatter¬ 
ing  of  molecular  hydrogen  from  Ag(lll)  is  a  prime  example  for 
which  to  attempt  an  empirical  analysis  of  the  physisorptive  inter¬ 
action  potential.  The  extremely  low  corrugation  of  this  surface 
facilitates  the  analysis  of  DSA  resonances  since  complications 
due  to  band  structure  are  not  present.  The  smoothness  of  this 
surface  also  concentrates  the  final  scattering  intensities  in 
the  specular  channel.  Molecular  hydrogen  does  not  dissociatively 
chemisorb  on  Ag(lll) ,  thereby  allowing  cryogenically  cooled  targets 
to  be  used  and  minimizing  Debye-Waller  attenuation.  H2,  D2 , 
and  HD  are  excellent  candidates  for  study  in  the  sense  that  only 
a  small  number  of  final  inelastic  scattering  channels  are  energe¬ 
tically  accessible.  Theoretically,  this  also  makes  potential  vari¬ 
ation  within  the  context  of  close  coupled  calculations  feasible, 
which  enables  one  to  take  full  advantage  of  the  information  in  the 
selective  adsorption  resonances  and  inelastic  transition  probabil¬ 
ities.  The  relatively  large  rotational  spacings  of  these  molecules 
also  allow  the  final  rotational  state  populations  to  be  probed  by 
measuring  outgoing  angular  distributions,  alleviating  the  need  for 
a  laser-based  detection  system.  Finally,  it  is  also  worth  noting 
that  although  HD/Ag(lll)  is  the  richest  system  in  terms  of  large 
rotational  inelasticity  and  producing  strong  RMSA  resonances,  H2 
and  D2  actually  provide  a  more  sensitive  probe  of  the  anisotTopy 
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of  Vqq(z,  0).  This  interesting  situation  arises  since  the  HD 
rotational  transitions  and  J-dependent  RMSA  energy  shifts  are 
primarily  governed  by  the  mass  asymmetry  of  the  molecule,  which 
is  large  enough  to  obscure  the  subtle  effects  which  are  due  to 
the  weak  anisotropy  of  the  potential.  Such  problems  do  not  arise 
when  similar  measurements  are  performed  with  the  homonuclear  iso¬ 
topes  H2  and  D2- 

The  problem  of  inversion  of  a  multiple  component  potential 
such  as  V0Q(z,  0)  presents  a  major  theoretical  challenge.  In  the 
absence  of  good  inversion  procedures  we  adopt  here  an  empirical 
approach,  taking  a  suitable  parametric  form  for  V0Q(z,  e)  and  then 
performing  quantum  close  coupled  scattering  calculations  to  obtain 
transition  probabilities  and  resonance  shifts  in  agreement  with 
experiment.  The  rotational  transition  probabilities  and  the  reson¬ 
ance  energies  provide  complementary  information  on  the  interaction 
potential  which  will  prove  to  be  essential  for  an  unambiguous  deter¬ 
mination  of  the  anisotropy  of  the  interaction.  We  proceed  stepwise 
by  assuming  the  previously  established  isotropic  component  Vqq(z) 

[6,  10]  and  consider  here  only  the  remaining  parameterization  of 
the  anisotropic  component. 

Finally,  hydrogenic  systems  are  also  simple  enough  that 
ab  initio  theoretical  estimates  of  the  interaction  potential  can 
be  made  and  asymptotic  expressions  for  both  attractive  and  repulsive 
terms  of  the  anisotropic  component  of  Vqq(z,  0)  are  available  for 
comparison  [12,  13]. 

The  plan  of  the  remainder  of  the  paper  is  as  follows.  In 
Section  II  we  discuss  the  theoretical  basis  for  a  single  functional 
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form  of  the  laterally  averaged  potential  and  outline  the  theo¬ 
retical  procedure  to  be  followed  in  the  investigation  of  the 
potential.  Section  III  contains  a  summary  of  the  experimental 
procedures  for  obtaining  rotationally  inelastic  transition  prob¬ 
abilities,  resonance  energies  and  associated  linewidth  informa¬ 
tion.  In  Section  IV  the  rotationally  inelastic  scattering  data 
and  data  analysis  for  ^  and  are  presented,  followed  by  the 
results  of  scattering  calculations.  In  Section  V  the  selective 
adsorption  resonances  for  n-I^,  P-^,  n-D2  and  o-D^  are  discussed: 
linewidth  broadening  of  n-H2  due  to  averaging  over  |J,m>  components 
of  the  initial  beam  and  shifts  of  n-H2  resonances  relative  to  p-H2 
are  analyzed.  The  analysis  in  each  of  these  sections  yields  a 
set  of  possible  parameteri zations  of  the  anisotropy  and  the  inter¬ 
section  of  these  two  sets  gives  us  a  unique  parameter  set.  In 
Section  VI  we  present  the  rotational  transition  probabilities  for 
the  HD  scattering  and  compare  the  HD  RMSA  resonance  energies  and 
transition  probabilities  with  the  theoretical  predictions  from 
the  best  fit  anisotropy  obtained  from  the  H2  and  D^  data.  Finally, 
Section  VII  contains  a  discussion  of  the  results,  a  summary,  and 
our  concluding  remarks. 

% 

II .  THEORY 

A.  The  Laterally  Averaged  Interaction  Potential 

The  laterally  averaged  interaction  potential  for  a  rigid 
rotor,  VQ0(z,  0),  may  be  written  as  the  sum  of  an  attractive  long 
range  polarization  potential  and  a  repulsive  short  range  potential 
deriving  from  overlap  of  molecular  and  solid  electronic  wave- 
functions.  In  general,  each  of  these  terms  will  have  a  different 
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orientation  dependence: 

v00(z-  e>  ■  V00,r'z-  9)  *  V00,A(z-  9)  <2-» 

Ab  initio  theoretical  calculation  of  each  of  these  components 
is  rendered  exceedingly  difficult  by  the  many  body  nature  of  the 
interaction  and,  for  metals,  also  by  the  delocalization  of  the 
electron  density.  However,  asymptotic  expressions  have  recently 
been  derived  which  givi.  rise  to  simple  few  parameter  functional 
forms.  For  the  van  der  Waals  interaction 

-C,(6) 

V00,A^Z>  9^  *  T  '2'2> 

where  the  orientation  dependence  of  0^(0)  arises  as  a  result  of 
the  anisotropy  of  the  molecular  polarizability.  For  a  homonuclear 
diatomic  molecule  it  is  given  by  [12] 

C3(6)  =  cj°>  +  C^2)  P2(cose)  (.2.3) 

Harris  and  Feibelman  have  obtained  C^2^  ^  0.05  for  H2  [12]. 

The  repulsive  component  has  been  estimated  by  Harris  and  Liebsch 
from  consideration  of  the  effective  orthogonalization  of  the  mole¬ 
cular  and  metal  electron  wavefunctions  as  [13] 

-k  z  r2  k2 

V00,R(2’  e)  *  Ae  *  C1  +  “f T1  (1  +  2  p2tcose^  (2.4) 

2  2 

where  k$  =  (2  )  is  the  decay  constant  of  a  Fermi-surface 

Bloch  wavefunction  ($  is  the  workfunction,  mg  the  electron  mass) 

and  r  the  molecular  internuclear  distance, 
e 

These  estimates  give  a  laterally  averaged  interaction  poten¬ 
tial  with  two  parameters  8^  and  8^  specifying  the  anisotropy 


V00tz>  9)  ■  v00,Rtz)  [1  *  SR  t,2tc°se’]  *  V00,A(Z) 

[i  +  eA  P2(cose) ] 

=  vSo(z)  +  voo(z:>  p2Ccos0) 


(2.5) 

(2.6) 


Using  the  appropriate  parameters  for  H^/Agflll)  (rg  =  1.4  a^  , 

$  =  4.24  eV  we  find  0R  ^  0.10  from  Equation  (2.4)  and  0A  %  0.06 
from  equation  (13)  of  ref.  [12],  using  accurate  values  of 
the  polarizability  tensor  of  H2  [14].  Both  0A  and  0R  are  posi¬ 
tive,  reflecting  the  asymptotically  preferred  perpendicular  orien¬ 
tation  at  large  z  and  parallel  orientation  at  small  z.  Equations 
(2.5)  and  (2.6)  form  the  basis  of  the  parameterizations  used  in 
our  studies.  We  shall  use  two  different  isotropic  potentials 


v00(-2-)  v00 ,  Z^  +  V00,A^Z)  (2,7) 

2 

and  obtain  the  anisotropic  term  vQ0(z)  from  Equation  (2.5)  as  the 
linear  combination 


V00^Z^  6R  V00,R(‘Z^  +  SA  V00,A^  t2,8) 

Thus  0A  and  0R  can  be  regarded  as  parameters  which  allow  independ¬ 
ent  variation  of  the  attractive  and  repulsive  components  respect¬ 
ively,  with  regard  to  the  molecular  orientation.  Equation  (2.8)  is 
an  improvement  over  the  one  parameter  anisotropy  previously  used  to 
study  rotationally  inelastic  scattering  of  HD  from  Ag(lll)  [9].  Th 
anisotropy  is  thus  determined  hoth  by  the  parameters  0A,  0R  and  by 
the  functional  form  of  the  isotropic  potential  Vqq(z) . 

For  Vqq(z)  we  use  (i)  the  variable  exponent  potential  (here¬ 
after  VEP) 


vSn  (2)  •  D  [(1  +  £  (z-z  J)'2p  -2(1  +  4  (z-zJ)’P] 


(2.9) 
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and  (ii)  the  exponential-3  potential  (hereafter  exp- 3) 


Daz 


-a(z-zj 


00 


,  ,  "WAe  ,r  3  eJ  ‘e.3 

(*)  =  -  f~ >  fC2^ 


f(z)  =  1-  [2(z/a0) (z/aQ  +  1)  +  1]  exp(-2z/aQ) 


(2.10) 

C2.ll) 


where  f(z)  is  a  cut-off  function  going  smoothly  from  f(0)  =  0 
to  f(z)  *  1  at  large  z.  The  VEP  is  a  flexible  three  parameter 
potential  with  semi-analytic  eigenvalues  [15]  which  has  been  used 
to  fit  the  bound  state  levels  of  laterally  averaged  physisorption 
potentials  for  atom-surface  systems  [16].  The  exponential-3  poten¬ 
tial  has  the  form  predicted  by  the  ab  initio  theoretical  studies, 
Equations  (2.2-4),  but  is  somewhat  more  difficult  to  work  with, 
having  no  analytic  eigenvalues  and  requiring  a  cut-off  function  to 
remove  the  singularity  at  z  =  0 . 

Equation  (2.7)  with  both  (2.9)  and  (2.10)  is  approximate  in 

that  it  neglects  correlation  between  metal  and  electronic  wave- 

functions.  A  consequence  of  the  use  of  the  asymptotic  results 

for  Vqq.  a(z,  6)  and  VQ0  R(z,  6)  (Equations  (2.2-4))  is  that  the 
*  2 

z  dependence  of  Vqq(z)  is  determined  by  that  of  the  components 
v00  R(z)  and  vQ()  A(z)  of  Vg0(z).  The  calculations  presented  in 
Section  IV  will  show  that  this  form  cannot  entirely  account  for 
the  energy  dependence  of  the  transition  probabilities.  An  analog¬ 
ous  situation  is  seen  in  diffractive  scattering  of  atoms  from 
surfaces  [17],  We  have  therefore  also  briefly  considered  the 
introduction  of  a  third  parameter  into  Vqq(z)  which  modifies  the 
softness  of  the  repulsive  component: 


(z) 


"  SR  e*P(~a2z)  v00,Rfz)  +  6  A  V00,A(z) 


(2.12) 


BmijsnrarauiMAMflM  man jjjj m xB&xxt,  a m>w^  J5H 
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The  above  discussion  based  on  Equation  (2.5)  refers  to 
the  interaction  potential  for  a  homonuclear  molecule.  The  inter¬ 
action  potential  for  HD  in  the  center  of  mass  coordinates  (z,  6) 
is  then  obtained  from  that  for  H2  (Equation  (2.5))  by  the  trans¬ 
formation 

z  =  z'  -  6cos6  (2.13) 

with  Equation  (2.5)  valid  for  (z*,  6).  Thus  for  HD  we  obtain 

VQ0(z,  6)  *  VQ0  ( z  * ,  0)  (2.14) 

*  v®0(z  +  <5  cos 6)  +  Vq q ( z  +  6  cose)  P2(cos0)  (2.15) 

v00(z)  +  l  -  0,1^2,...  v00(z)  pACcos0)  (2.16) 


with  6  ■  r  /6,  and  where  r  is  the  equilibrium  internuclear  dis- 

6  v 


tance.  The  resulting  potential,  denoted  VQ0(z,  0)  to  distinguish 
it  from  the  homonuclear  case,  now  contains  odd  oTder  Legendre 
terms,  reflecting  the  absence  of  the  AJ  *  0,2  selection  rule  for 


the  asymmetric  species.  The  separation  of  the  zeroth  order  term 


into  vjjg(z)  and  vjjg(z)  is  convenient  in  order  to  be  able  to  refer 


to  the  resonance  states  (both  DSA  and  RMSA)  of  all  isotopic  vari 
ants  in  the  same  terms  (Sections  V  and  VI) . 


B.  Theoretical  Analysis  of  Transition  Probabilities  and  Resonance 
Uatl - -  - 


It  is  clear  that  rotational  transition  probabilities  and  selec¬ 


tive  adsorption  resonance  data  provide  independent,  and  in  some  sense 

.2 


complementary,  information  on  the  anisotropic  component  Vq0(.z)  of 

the  interaction  potential.  The  rotational  transition  probabilities 

probe  the  anisotropic  component  at  small  z  values  near  the  turning 

.0 


point  of  Vqq(z)  and  will  be  determined  both  by  BR  and  B^-  At  high 
incident  energies  the  BR  dependence  may  be  expected  to  be  stronger 


than  the  8A  dependence  but  in  general  the  precise  balance  will 
be  dependent  on  the  form  of  the  isotropic  potential  components 
Vqq  r(z)  and  a( z ) .  The  resonances  are  determined  primarily 

by  the  isotropic  potential  Vqq(z):  information  on  the  anisotro- 
pic  component  in  the  region  of  the  bound  levels  of  Vq0(z)  is 
obtained  from  the  |Jm>  -dependent  shifts  and  linewidth  differences 
seen,  for  example,  for  n-H2  and  p-H2.  Since  the  form  of  vttQ(z) 
in  the  region  of  the  minimum  of  VqqCz)  strongly  dependent  upon 
the  combination  of  8A,  6R*  these  shifts  and  linewidths  will  be 
very  sensitive  to  simultaneous  variations  in  both  eA  and  8R. 

We  first  discuss  the  rotational  transition  probabilities. 

The  rotational  inelasticity  of  HI)  is  considerably  larger  than 

that  of  H2  and  D2  as  a  result  of  the  asymmetry  introduced  by  the 

offset  of  the  center  of  mass,  5  (Equation  (.2.13)).  However,  as 

a  consequence  the  dependence  of  rotational  inelastic  probabilities 

on  the  relatively  weak  orientational  anisotropy  deriving  from 
2 

Vqq(z)  is  masked  by  the  effect  of  this  center  of  mass  asymmetry. 

Far  more  sensitive  to  BA  and  (JR  are  the  H2  and  D2  rotational  tran¬ 
sition  probabilities  as  a  function  of  incident  perpendicular  energy 
2 

Ez  *  Ejj  cos  8^  where  E^  is  the  heam  energy  and  0^  is  the  incident 
beam  angle  from  the  surface  normal.  The  experimentally  measured 
ratio  of  probabilities  P(J  =  0  -*•  2)/P(J  ■  0  -*■  0)  for  H2  and  D2  is 
thus  the  quantity  which  will  be  compared  with  the  transition  prob¬ 
ability  ratio  derived  from  close  coupling  calculations  in  ordeT  to 
extract  a  first  quantitative  estimate  of  8A>  6R  for  a  given  v^Q(z). 
Note  that  the  relatively  large  HD  rotational  transition  probabili¬ 
ties  are  nevertheless  dependent  upon  gross  features  of  the  isotro¬ 
pic  potential  vj^fz),  such  as  the  well  depth.  This  was  utilized 


in  our  previous  study  of  the  isotropic  component  of  Vqq(z,  e) 

do]. 

In  the  first  approximation  DSA  resonance  energies  are  given 
by  eigenvalues  of  Vqq(z)  plus  the  kinetic  energy  parallel  to  the 
surface  obtained  by  transfer  of  a  reciprocal  lattice  vector  G  to 
the  incident  parallel  momentum  K^.  Similarly,  RMSA  resonance 
energies  are  given  to  a  first  approximation  by  the  bound  vibra¬ 
tion-rigid  rotor  (BVRR)  energies  [18],  which  are  eigenvalues  of 
VqqCz)  plus  free  rotor  energies.  Shifts  from  these  zeroth  order 
DSA  resonance  energies,  and  from  the  BVRR  energies,  arise  from 
the  anisotropic  component  of  Vqq(z,  0)  and  result  in  the  experi¬ 
mentally  observed  jJm>  and  J  dependence  of  bound  levels  extracted 
from  DSA  and  RMSA  measurements,  respectively. 

The  connection  between  the  two  kinds  of  resonances  is  clearly 
seen  from  the  standard  set  of  close  coupled  equations  resulting 
from  a  product  basis  expansion  in  rigid  rotor  and  diffraction 
states  with 


X1  2 


H  =  TR  v  “  zr  [iTne  36  sin0  30]  +  VQ0Cz>  0) 


(2.17) 


<Kz,0,<j>,R)  *  l  X(z)  YT  (0,0)  e 
GJm  GJm  Jra 


i(K.  +  G) »R 


(2.18) 


we  obtain 


X(z)  =  I  B  <Jm|VQO(z,0)|J’m>  x  (z) 
GJm  J'dT  uu  GJ'm 

k2GJ  “Jr  E  -  +  9z  -T  JCJ  +  « 

~  'K 


(2.19) 


(2.20) 


For  atoms  or  spherical  molecules  Vqq(z,0)  reduces  to  VqQ(z)  and 
the  resonance  condition  is  then 
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2  2 

E  -  (*i  +  9 2  +  TT  JCJ  +  +  en  £2.21) 

with  en  a  negative  eigenvalue  of  Vqq(z).  Equation  C2.21)  is  the 
generalization  of  the  zeroth  order  atom-surface  DSA  resonance 
condition  and  can  be  regarded  as  representing  a  diffractive  and 
rotationally  excited  molecular  state  in  which  the  molecule  moves 
freely  parallel  to  the  surface  but  is  in  a  bound  vibrational 
state  perpendicular  to  the  surface  with  energy  en  <  0 .  We  refer 
to  these  energies  as  BVRR  energies,  for  both  zero  and  nonzero  G. 

In  general  the  set  of  coupled  Equation’s  C2.19)  must  be 
solved  exactly  for  each  G  and  m  in  order  to  find  the  resonance 
energies.  However,  when  the  anisotropy  is  small  the  resonances 
may  be  found  hy  perturbation  theory  [18].  Thus  for  the  homo- 
nuclear  diatom-surface  systems,  application  of  first  order  per¬ 
turbation  theory  to  Equation  (.2.19)  gives  resonance  energies 


E  =  (*i 


G)2 


TT  J(J  + 


'n 


Cl) 

'n,  Jm 


where 


(2.22) 


,(D 


<  n 


(2.23) 


while  for  the  heteronuclear  case  we  have  Equation  (2.22)  with 

rePlaced  by 

*  l  <n|v?0(z)  |n><Jm|PpCcos6)  | Jm>  (.2.24) 

n’Jm  l  -  0,1,2... 00  1 

Equation  (2.22)  is  the  kinematic  resonance  condition  for  a  mixed 
DSA-RMSA  resonance  to  first  order  in  the  anisotropic  component  of 
V00(z»  e)*  Pure  RMSA  or  pure  DSA  resonances  are  obtained  by  set¬ 
ting  G  or  J  equal  to  zero  respectively  (if  the  initial  state  is 


rotationally  excited,  pure  DSA  is  obtained  when  AJ  =  0) . 

^■n^Jm  an<*  ^n^Jm  are  t^ie  ^rst  order  |Jm>  -  and  n-dependent 

shifts  from  the  BVRR  energies  (Equation  (2.21))  for  the  homo- 

nuclear  and  heteronuclear  adsorbed  molecular  species  respectively. 

Equations  (2.22)  and  (2.24)  show  that  even  in  the  absence  of  an 

2 

anisotropic  term  vqq(z),  the  observed  HD  RMSA  resonance  energies 
will  show  |Jm>  -dependent  shifts  for  a  given  bound  state  n.  The 
fact  that  the  observed  J-dependent  shifts  seen  in  HD  RMSA  reson¬ 
ance  measurements  on  both  Pt(lll)  [8]  and  Ag(lll)  [9]  (for  m  *  0) 
are  much  larger  than  the  |Jm>  -dependent  shifts  seen  in  H2  and  D2 
DSA  and  RMSA  measurements  on  Ag(lll)  [6]  (^l.OmeV  for  HD,  compared 
to  ^<0.2meV  for  H2  and  D2)  is  not  surprising  in  view  of  the  large 
asymmetry  of  HD.  Thus  just  as  with  the  rotational  transition  prob¬ 
abilities,  the  high  ailisotropy  induced  by  the  asymmetry  of  HD  masks 

2 

the  potential  anisotropy  from  Vqq(z)  and  renders  the  RMSA  resonance 

2 

shifts  for  HD  less  sensitive  to  small  changes  in  Vq0(z)  than  the 
DSA  and  RMSA  resonance  shifts  for  H2  and  D2- 

Provided  first  order  perturbation  theory  is  a  valid  descrip¬ 
tion  of  the  resonances,  then  an  adsorbed  ground  state  molecule 
|Jm>  *  1 0 0 >  will  have  zero  shift  (or  S^Jm^  and  yields  the 

true  bound  state  energies,  en,  of  VqQ(z) .  This  forms  the  basis  of 
the  empirical  determination  of  the  isotropic  component  of  VQ0(z,  0) 
from  DSA  measuremer :s  of  p-H2  and  o-D2  on  Ag(lll)  which  was  reported 
previously  [10  J.  (Note  that  adsorbed  J  *  0  can  not  be  observed  in 
RMSA  resonances,  by  definition).  We  should  emphasize  that  because  the 
bound  state  energies  determined  for  |J,m>  *  J 0 , 0>  are  independent  of 
the  anisotropic  potential,  Vqq(z)  may  be  determined  first,  and  then 
the  anisotropic  components  may  be  determined,  without  the  need  for 
iteration. 


?  k^i  '  '  m  ’  if  ^  j  Fi  nar  -jr  ,’  _* -jt r > 

For  a  given  bound  level  n,  the  first  order  resonance  shifts 

en^Jm  Prov^^e  direct  information  on  the  anisotropic  component, 

2  * 

Vqq(z)  .  Any  well  resolved  selective  adsorption  resonance  which  is 


unambiguously  assignable  to  a  single  |Jm>  rotor  state  will  give 
rise  to  a  set  of  possible  pairs  of  parameters  BA,  BR  via 
Equations  (2.8)  and  (2.23).  Resonances  observed  for  n-H2  and 
n-D2  will  reflect  an  average  over  |Jm>  states  |0,0>,  |1,0>  and 
|l,±l>  (±  m  have  equal  shifts  according  to  equation  (2.23)), 
with  appropriate  weightings  resulting  from  the  nuclear  statistics. 
Unfortunately  the  resolution  in  the  present  experiments  (^25QyeV) 
is  too  low  to  resolve  the  individual  multiplet  components  for  the 
normal  hydrogen  isotopes  in  the  extremely  weak  DSA  resonance  dips 
observed  in  the  specular  scattering  from  Ag(lll) .  On  the  more 
strongly  corrugated  Ag(110)  surface  where  the  resonance  structure 
is  stronger,  there  is  some  evidence  for  a  splitting  of  the  n-H2 
DSA  resonances  [7,  19].  In  our  experiments  on  Ag(lll)  we  do,  how¬ 
ever,  observe  (i)  an  averaged  effective  shift  of  the  resonance  energy 
for  n-H2  from  that  for  p-H2  (and  n-D2  from  o-D2)  and  (ii)  a  dis¬ 
tinct  linewidth  broadening  for  n-H2  and  n-D2  relative  to  p-H2  and 
o-D2  respectively  [6].  Both  the  effective  shift  and  the  linewidth 
broadening  may  be  simulated  as  a  function  of  the  parameters  8A»  BR, 
using  a  natural  linewidth  for  the  |Jm>  states  obtained  from  the 
measured  linewidth  of  p-H2  and  o-D2» 

Each  of  these  comparisons  between  experiment  and  theoretical 
simulation,  i.e.  for  transition  probabilities  and  for  resonance 
shifts,  will  provide  a  set  of  parameter  pairs  6A»  SR.  We  regard 

the  intersection  of  these  two  sets  as  providing  the  correct  para- 

2 

meters  for  the  anisotropy  Vqq(z),  derived  from  a  given  isotropic 
potential  Vqq(z)  according  to  Equations  (.2.5)  and  (.2.6).  Because 
of  the  smaller  sensitivity  of  both  the  transition  probabilities 
and  RMSA  resonance  energies  for  HD  to  the  form  of  VqqCz) ,  we  shall 


use  this  as  a  check  to  differentiate  between  extremely  different 

regions  of  the  0A,  0R  parameter  space  but  not  to  investigate  small 

changes  in  the  values  of  8^  and  0R.  This  is  an  improvement  on  our 

earlier  study  of  the  HD/Ag(lll)  system  [9]  where  for  a  much  cruder 
2 

model  of  Vqq(z),  namely  0  =  0A  *  0R,  analysis  of  the  RMSA  resonance 
energies  alone  gave  a  preliminary  order  of  magnitude  estimate  of  0 
(representing  an  average  of  0A  and  0R) . 

C .  Close  Coupled  Calculations 

The  close  coupled  equations  to  be  solved  are  Equations  C2.19) 
for  G  =  0  and  m  =  0  (initial  rotor  state  |Jm>  ■  |00>).  The  total 

<v  ' 

energy  E  is  then  replaced  by  Ez  ,  the  z  component  of  the  incident 
beam  energy.  For  HD,  VQ0(z,  8)  is  replaced  by  VQ0(z,  8),  Equations 
(2.14-16).  The  use  of  the  flat  surface  assumption  to  obtain  the 
rotationally  inelastic  transition  probabilities  is  justified  for 
HD  by  the  extremely  small  size  of  the  diffraction  relative  to  the 
rotational  excitation,  and  for  H2  and  D2  by  the  absence  of  diffrac¬ 
tion  peaks  for  measurements  made  along  the  azimuth  of  weakest  cor¬ 
rugation.  To  ensure  an  accurate  comparison  of  resonance  positions 
with  the  experimental  data  we  replace  the  rigid  rotor  energies  by 
the  empirical  rotor  energies  E(J)  (Table  1) .  For  HD  at  beam 
energies  E^  'v*  120meV  we  have  only  four  open  rotational  channels. 

For  H2  and  D2,  we  have  four  and  six  open  channels,  respectively, 
which  are  reduced  by  symmetry  to  two  open  channels  for  H2  scat¬ 
tering  from  initial  state  J  *  0  and  3  open  channels  for  D2 .  The 
close  coupled  equations  are  solved  with  the  variable  step-variable 
interval  size  (VIVAS)  integrator  [20],  Typical  running  times  are 
6  sec.  C.P.U.  time  for  a  four  channel  system  at  E  *  60meV  on  a 


Vax  11/780  computer.  The  speed  of  this  integrator  is  extremely 
useful  for  the  multiple  runs  needed  for  an  iterative  search  of 
the  8^,  gR  parameter  space  even  though  the  further  time  saving 
usually  available  in  multiple  calculations  at  different  energies 
is  not  possible  when  potential  parameters  are  varied. 

The  best  fit  parameters  for  the  VEP  and  exp- 3  forms  of 
Vqq(z),  Equations  (2.9)  and  (2.10)  are  given  in  Table  2  [6,10]. 

The  first  order  resonance  shifts  for  H~  and  D0  were 

n,  Jin  2  2 

evaluated  directly  from  Equation  (2.23),  using  a  Numerov-Cooley 
algorithm  to  evaluate  the  bound  state  eigenfunction  with  energy 
en  and  then  integrating  over  vQ0  R(.z)  and  vQ0  ACz)  (Equation 
(2.8))  with  a  6300  point  Simpson’s  rule  quadrature. 

III.  EXPERIMENTAL  PROCEDURE 

The  experiments  described  in  this  paper  were  carried  out  in 
an  ultra-high  vacuum  beam-surface  scattering  apparatus  which  has 
been  previously  described  [9,21].  It  consists  of  the  following 
sections:  a  differentially  pumped  beam  source  manifold,  a  UHV 
crystal  chamber  which  includes  an  Auger  spectrometer  and  an  ion 
sputtering  gun,  a  douhly-dif ferentially  pumped  rotatable  mass 
spectrometer,  and  a  PDP-11  minicomputer  which  controls  the  experi¬ 
mental  apparatus  and  collects  the  data  through  a  CAMAC  interface. 

Tire  P-H2  and  0-D2,  i.e.  J  ■  0  hydrogen,  used  in  these  studies 
were  generated  by  passing  either  n-l^  or  n-D2  through  an  in-line 
catalytic  converter.  This  device  essentially  consists  of  a  co¬ 
axial  stainless  steel  tube  assembly  which  contains  a  dispersed 
nickel  catalyst  in  its  outer  section.  During  operation  the  con¬ 
verter  is  immersed  in  a  liquid  hydrogen  dewar.  This  induces  the 
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nuclear  spin  state  distribution  to  relax  to  its  equilibrium  value 
at  20  K,  i.e.  99.81  P*^  and  97.8%  A  detailed  description 

of  this  converter  has  been  previously  presented  in  our  earlier 
paper  dealing  with  the  hydrogen/Ag( 111)  isotropic  potential  [ 10 3 . 
The  beam  parameters  used  in  this  study  were  as  follows:  incident 
beam  divergence  0.1°,  Av/v  *  4.5%  (H2)  and  6.5%  (D2) • 

The  HD  used  in  these  experiments  was  also  generated  cataly- 
tically  by  passing  an  8:1  H2^D2  mixture  through  a  stainless  steel 
column  which  was  packed  with  magnesium  granules  and  heated  to  an 
operating  temperature  of  600  K.  This  8:1  feed  gas  ratio  was 
chosen  in  order  to  produce  an  HD  beam  of  high  intensity  and  narrow 
velocity  spread  (Av/v  ^  5%)  [8,9,22].  An  in-line  liquid  nitrogen 
cooled  molecular  sieve  trap  was  placed  between  the  HD  generator 
and  the  beam  source  in  order  to  remove  any  trace  moisture  from 
the  resulting  4:1  H2:HD  gas  mixture.  The  HD  after  supersonic 
expansion  is  predominantly  (>  991)  in  the  J  -  0  state.  This  is 
implied  from  the  absence  of  energy  gain  peaks,  J  =  1  0 ,  in  the 

experimentally  observed  angular  distributions  of  HD  scattering 
from  the  silver  surface.  (When  the  quality  of  the  jet  expansion 
is  intentionally  degraded  by  lowering  the  stagnation  pressure 
J  -  1  -*•  0  transitions  can  be  clearly  seen.)  The  HD  is  typically 
expanded  through  a  20  micron  nozzle  at  a  pressure  of  260  PSI. 

Three  types  of  scattering  data  are  reported  in  this  paper: 
angular  distributions  which  include  rotationally  inelastic  peaks, 
time-of- flight  distributions,  and  selective  adsorption  scans.  The 
angular  distributions  were  obtained  by  scanning  the  detector  at 
0.2°  increments  while  holding  the  angle  of  incidence  at  a  fixed 
value.  During  these  measurements  the  incident  beam  was  modulated 


with  a  50%  duty  cycle  chopper.  This  enabled  two  32-bit  scalers 
which  were  synchronized  to  the  chopper  to  collect  signal  plus 
background  and  background  only  information.  The  appropriate  sub¬ 
tractions  were  made  after  six  seconds  of  counting  at  each  0.2° 
angular  increment.  Time-of- flight  measurements  were  obtained 
with  a  1%  duty  cycle  dynamically  balanced  8-slot  chopper  (0.020 
inch  wide  slots,  6  inch  diameter,  typically  running  speed  of 
200-350  Hz)  operating  in  conjunction  with  a  255  channel  multi¬ 
channel  scaler  (0 . 25ysec/channel  dwell  time).  Selective  adsorp¬ 
tion  data  was  collected  by  measuring  the  peak  intensity  of  a 
given  rotational  and/or  diffractive  peak  at  incident  angle  incre¬ 
ments  of  0.2°.  In  practice,  these  data  were  collected  by  alter¬ 
nately  moving  the  crystal  polar  angle  by  ca.  0.2°  and  then  scan¬ 
ning  the  detector  through  the  desired  peak  with  the  computer 
recording  in  real  time  both  the  intensity  and  corresponding  in¬ 
cident  and  reflected  polar  angles.  Peak  heights  determined  in 
this  manner  include  contributions  from  both  the  coherently  scat¬ 
tered  molecules  and  the  diffuse  inelastic  background. 

The  Ag(lll)  crystal  was  cut  from  a  99.999%  boule  and  initially 
polished  at  the  Cornell  Materials  Preparation  Laboratory.  It  was 
subsequently  repolished,  aligned  with  Laue  X-ray  backreflection , 
and  mounted  using  procedures  that  were  previously  described  in 
detail  [9].  Crystal  quality  was  confirmed  by  observing  uniform 
and  high  He  and  H2  reflectivities  at  several  points  along  the 
surface  of  the  crystal. 


IV.  h2  and  d2  rotationally  inelastic  scattering 

A.  Experimental  Results 

In  this  section  the  experimentally  determined  rotationally 
inelastic  J  =  0  -*■  2  transition  probabilities  are  presented  for 
both  H2  and  D2*  The  convention  used  for  labelling  the  azimuthal 
angle,  <p ,  is  to  assign  <j>  =  0°  to  the  <112>  direction  and  <j>  =  30° 
to  <0ll>.  The  rotationally  inelastic  data  presented  in  this  sec¬ 
tion  were  taken  at  <J>  *  15°  to  avoid  the  overlaps  that  occasionally 
occur  in  the  final  angular  distributions  between  weak  diffraction 
peaks  and  weak  J  =  0  -*■  2  transitions.  This  overlap  of  diffractive 
and  inelastic  transitions  was  most  severe  along  =  0°.  As  a 
check  of  our  procedure,  inelastic  peaks  that  could  be  clearly 
measured  at  <j>  =-  0  were  fully  analyzed,  and  found  to  be  in  agree¬ 
ment  with  the  (j>  =  15°  data  within  experimental  uncertainty.  The 
diffraction  and  simultaneous  rotation-diffraction  intensities 


will  he  discussed  in  a  later  publication  [11]. 

Figure  2  shows  the  angular  distributions  for  n-H2  scatter¬ 
ing  from  Ag(lll)  at  four  incident  polar  angles  which  are  separated 
sequentially  by  5°  angular  increments.  The  J  =  0  -*■  2  peak  is 
barely  observable  at  9^  *  40°  and  gradually  emerges  as  a  distinct 
peak  in  the  shoulder  of  the  J  *  0  -*■  0  specular  peak  as  0^  decreases 
towards  the  surface  normal  (i.e.,  the  inelastic  transition  prob¬ 
ability  increases  relative  to  the  specular  beam  as  E  is  system- 

zi 

atically  increased) .  The  FWHM  of  the  specular  beam  angular  pro¬ 
files  shown  in  Fig.  2,  vl.0°,  are  just  somewhat  larger  than  that 
of  the  incident  beam,  ^0*7°,  as  expected  for  a  high  quality  sur¬ 
face.  The  FWHM  of  the  J  ■  0  -*■  2  peaks  are  of  course  larger  due 
to  the  kinematics  of  the  scattering  process  -  different  velocities 


present  in  the  incident  beam  velocity  distribution  are  inelas- 
tically  scattered  into  different  outgoing  angles. 

Figure  3  shows  the  angular  distributions  for  n-D2  scatter¬ 
ing  from  Ag(lll)  for  a  similar  sequence  of  four  incident  polar 
angles  (out  of  a  total  of  six  angles  measured  for  n-D^) .  The 
same  trend  of  increasing  inelastic  transition  probabilities  as 

0.  decreases  towards  normal  incidence  (E  increasing)  is  again 
x  2i 

present.  Also  note  that  the  inelastic  peaks  for  D2  are  scaled 
up  by  a  factor  of  10  relative  to  the  specular  channel  while  those 
for  H2  are  scaled  by  a  factor  of  100  even  though  the  relative 
population  of  J  =  0  in  the  incident  beam  is  only  about  twice  that 
for  This  implies  that  the  relative  inelastic  transition 

probabilities  of  D2  are  significantly  larger  than  those  for  H 2» 
This  is  to  be  expected  from  the  much  smaller  rotational  energy 
spacing  of  relative  to  The  J  *  1  -*■  3  transition  is  also 

energetically  allowed  for  D2,  and  is  observed  at  several  angles 
(Figure  3) . 

B .  Experimental  Analysis 

The  final  scattering  angles  for  non- resonant  rotationally 
inelastic  transitions  are  given  by  the  following  two  kinematic 
constraints:  (i)  conservation  of  parallel  momentum 

Kf  -  Kf  +  G  (3.1) 

and  (ii)  conservation  of  total  energy 

k?  -  k^  +  2m  [E(jf)  -  E(J.)]  (3.2) 

■ff 


and  are  found  to  be  in  excellent  agreement  with  the  angles  cal¬ 
culated  using  empirical  rotor  energies  in  Equation  (3.2)  with 
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G  =  0 ,  i .e. 

k2  "  k*  +  %  -  E(J. )  3  (3.3) 

zi  zf  “fi  1  1 

This  confirms  that  only  the  perpendicular  component  of  the  beam 
energy  is  involved  in  the  energy  transfer  process.  This  is  demon¬ 
strated  in  Figure  4  which  shows  two  HD  angular  distributions  taken 
at  the  same  E2  but  using  considerably  different  incident  beam 
energies . 

In  order  to  directly  compare  the  experimental  results  with 
transition  probabilities  obtained  from  quantum  close  coupled  scat¬ 
tering  calculations  a  detailed  data  reduction  procedure  was  carried 
out.  Experimental  transition  probabilities  were  obtained  after 
(1)  deconvoluti'ng  the  instrumental  contributions  from  the  angular 
distributions,  (2)  correcting  for  Debye-Waller  attenuation,  and 
(3)  accounting  for  the  quantum  statistics  of  the  initial  rotational 
state  population  distributions. 

1.  Deconvolution  of  Final  Angular  Distributions 

The  supersonic  molecular  beam  distribution  is  modelled  by  a 
standard  shifted  Gaussian: 

p(v)dv  -  Nv3exp[-(v-v0) 2/a2]dv  (3.4) 

where  Vq  is  the  stream  velocity  and  a  the  Gaussian  width.  This 
can  be  readily  transformed  to  the  flight  time  distribution  in 
terms  of  number  density,  which  is  the  quantity  directly  deter¬ 
mined  by  our  detector.  The  measured  time-of-f light  data  of  the 
incident  beam  is  predominantly  a  convolution  of  the  true  particle 
flight  time  distribution  with  the  chopper  gating  function.  This 
gating  function  was  simulated  by  overlapping  the  beam’s  spatial 
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diameter  at  the  chopper  (0.5  mm)  with  a  moving  chopper  slit 
(width  0.51  mm).  Convolution  of  an  assumed  flight  time  distri¬ 
bution  with  the  gating  function  was  fitted  to  the  experimental 
distribution  (after  subtraction  of  the  ion  flight  time  in  the 
detector)  with  a  nonlinear  least  squares  fittings  routine  to 
obtain  the  beam  parameters. 

Angular  distributions  were  then  fit  with  a  simulation  pro¬ 
gram  based  on  the  fitted  incident  beam  characteristics,  scatter¬ 
ing  kinematics,  and  the  two-dimensional  detector  resolution  func¬ 
tion  which  was  determined  experimentally  with  the  incident  beam 
passing  through  a  0.05  mm  aperature  in  the  detector.  This 
Gaussian-shaped  function  has  an  in-plane  FWHM  of  0.66°  and  out- 
of-plane  FWHM  of  0.40°.  The  diffuse  background  was  modelled  with 
a  Gaussian  function  in  terms  of  (0r~6<O /0  .  where  6n  is  the 
Gaussian  origin  and  0w  the  width.  These  parameters  were  either 
assigned  or  freely  fitted  to  get  the  best  results.  A  nonlinear 
least  squares  fitting  routine  was  used  with  the  background  model 
to  fit  the  angular  distribution  data.  The  area  of  the  coherent 
part  of  the  peak  (i.e.  free  of  diffuse  background)  was  then  scaled 
to  the  area  of  the  simulated  peak  derived  from  the  incident  beam 
parameters,  scattering  kinematics,  and  the  detector  function. 
Inelastic  transition  probabilities  (still  needing  to  be  corrected 
for  Debye-Waller  attenuation  and  incident  ensemble  spin  statistics) 
suitable  for  comparison  with  scattering  calculations  were  obtained 
in  this  manner.  A  typical  fit  obtained  with  this  procedure  is 
shown  in  Figure  5. 
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2 .  Correction  for  Thermal  Attenuation 

The  surface  temperature  dependence  of  the  elastic  scatter¬ 
ing  of  H2,  *  an<*  and  also  the  rotationally  inelastic  scat¬ 

tering  of  HD  from  Ag(lll)  has  been  analyzed  in  detail  in  a  prev¬ 
ious  paper  [10].  The  Debye-Waller  attenuation  factor  for  a  ro¬ 
tationally  inelastic  peak  can  be  derived  from  the  simplest  treat¬ 
ment  for  atom-surface  scattering: 


*  eXpCi^  (Vos20i  *  D)Ts 


rE^os  ef  +  d 


E^cos  0^  +  D 


(3.5) 


where  0£  and  0£  are  related  by 

EicosZ0f  -  Eicos20i  +dS2[Ji(Ji  +  1)  -  Jf(J£  +  1)  ]/2m 


(3.6) 


Here  m  and  m  are  the  masses  of  the  incident  molecule  and  solid 
&  5 

atoms  respectively,  0D  is  tie  surface  Debye  temperature  and  D 
the  well  depth  of  the  laterally  averaged  interaction  potential. 

The  parameters  0^  and  D  extracted  from  temperature  dependent 
studies  are  given  in  Table  3,  where  it  is  apparent  that  a  signi¬ 
ficant  discrepancy  is  present  between  the  homonuclear  and  hetero- 
nuclear  results.  Other  evidence  for  the  well  depth  discussed  in 
our  previous  work  on  the  isotropic  potential  component  [10]  in¬ 
dicates  that  the  value  of  ,\/32meV  is  the  correct  one.  Since  the 
applicability  of  standard  Debye-Waller  theory  to  molecule-surface 
collisions  is  not  well  established,  the  meaning  of  this  difference 
hetween  the  isotopes  is  unclear.  We  therefore  regard  these  para¬ 
meters  as  an  empirical  means  of  correcting  the  experimental  tran¬ 
sition  probabilities  for  thermal  attenuation  due  to  the  motion  of 
the  target  surface. 
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3.  Initial  Rotational  State  Distributions  for  H2  and  D2 

The  rotational  state  distribution  of  the  incident  molecular 
beam  can  be  calculated  utilizing  the  energy  balance  equation  [23] 

!  kuV  w  *  7 mv2  * !  Vp  ♦  w  (3-7> 

where  TQ  is  the  nozzle  temperature,  Tp  the  parallel  temperature 
of  the  beam,  Tr  the  rotational  temperature,  m  the  molecular  mass, 
v  the  terminal  parallel  velocity,  and  Er(T)  the  average  rotational 
energy  at  temperature  T.  This  energy  balance  equation,  when  com¬ 
bined  with  the  assumption  that  the  ortho  and  para  spin  components 
of  the  initial  ensemble  do  not  interconvert  during  supersonic 
expansion,  allows  the  terminal  rotational  state  distributions  to 
be  obtained.  The  results  of  this  analysis  are  shown  in  Table  4 
for  n-H2,  P-H2,  and  n-D2. 

C .  Experimentally  Derived  Inelastic  Transition  Probabilities 

The  rotationally  elastic  and  inelastic  transition  probabilities 
for  H2  and  D2  are  given  in  Tables  5  and  6,  respectively.  The 
J  *  0  -*  0  transition  probabilities  were  calculated  by  combining 
the  specular  intensity  data  from  angular  distribution  measurements 
with  the  initial  population  distributions  presented  in  Table  4. 

It  was  assumed  that  the  elastic  scattering  probabilities  for  all  of 
the  initially  populated  rotational  states  were  equal.  This  is  a 
valid  assumption  for  the  hydrogen/Ag(lll)  system  at  the  collision 
energies  used  in  this  study.  As  a  final  check  of  our  data  analysis 


procedure  we  also  determined  the  P(J  *  0  ■+■  2)/P(J  *  0  -*■  0)  ratios 
for  H,  and  D-  at  selected  collision  energies  using  beams  of  pure 


p-H2  an<*  °"®2*  These  ratios  were  found  to  be  consistent  within 
experimental  error  with  the  values  obtained  using  n-H2  and  n-D2. 
The  P(J  =  0  2)/P(J  *  0  -►  0)  ratios  presented  in  Tables  5  and  6 
are  the  resulting  relative  transition  probabilities. 


D.  and  Rotational  Transition  Probabilities:  Comparison 

Between  Experimental  Results  and  Quantum  Scattering  Calculations 

We  now  compare  the  experimental  ratio  of  transition  prob¬ 
abilities  P(J  *  0  +  2)/P(J  »  0  0)  with  the  results  of  scatter¬ 

ing  calculations  for  both  VEP  and  exp- 3  based  forms  of  the  poten¬ 
tial  VQq(z,  0).  Figure  6  shows  the  calculated  ratio  Pc(J  •  0  -*•  2)/ 
PC(J  *  0  -*•  2)  as  a  function  of  $R  for  a  range  of  values  of  ba,  at 

a  constant  normal  incident  energy  E  ■  53.2meV.  The  isotropic 
Q  1  2 

potential  v0Q(z)  and  thus  the  hasis  for  vftft(z)  here  is  the  VEP. 

The  experimental  ratio  is  denoted  by  the  solid  horizontal  line  and 

the  error  bounds  indicated  by  the  dashed  lines.  The  ratio 

Pc(J  *  0  -*•  0)/Pc(J  =  0  -►  0)  shows  a  smooth  variation  as  a  function 

of  3a  and  of  8R:  we  note  that  for  certain  values  of  ba  there  are 

two  values  of  8R  giving  the  correct  transition  probability  ratio. 

Repeating  this  search  for  each  of  the  four  normal  incident  energies 

measured  for  H2  and  for  the  six  normal  incident  energies  for  D2 

(Tables  5  and  6)  produces  at  each  normal  incident  energy  E  a 
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whole  set  of  parameter  pairs  BA,  BR  which  will  fit  the  H2  or  D2 
transition  probability  ratios,  within  the  experimental  error 
bounds.  In  Table  7  we  show  the  combinations  SA,  BR  derived  by 
scanning  both  3A  and  BR  in  increments  of  0.01  for  the  VEP  para¬ 
meterization:  the  table  contains  the  value  of  8R  for  given  8A 
which  brings  the  probability  ratio  P(J  *  0  -*>  2)/PCJ  *  0  -*•  0) 
closest  to  the  experimental  values  for  each  energy.  Although  one 


or  more  values  of  BR  may  fall  between  the  experimental  error 
bounds,  for  this  grid  spacing  there  are  never  two  successive 
values,  i.e.  BR  =  B  and  BR  =  B  +  0.01  within  these  limits.  Where 
successive  values  of  8R  give  ratios  lying  above  and  below  the 
error  bounds,  an  interpolated  value  of  6  +  0.005  has  been  used. 

For  the  values  of  8A  showing  two  possible  values  of  BR  (e.g. 

Ba  =  ±  0.11  in  Figure  7)  simultaneous  comparison  at  different 
energies  allows  one  possibility  to  be  eliminated  because  of  an 
unrealistically  sharp  change  in  Bn  with  E„  .  Table  7  and  the 
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corresponding  results  for  the  exp- 3  parameterization  are  summarized 

in  Figures  7  and  8.  The  crosses  (X)  and  pluses  (+)  represent 

values  of  BA»  8R  which  fit  the  experimental  transition  probability 

for  H2  and  D2  respectively,  at  any  one  of  the  above  criterion. 

The  circles  in  Figure  8  represent  combinations  6A*  BR  derived 

from  the  n-H2  shifts  and  linewidth  broadening,  to  be  discussed  in 

Section  V.  Figures  7  and  8  thus  show  the  combinations  fiA»  BR 

which  are  most  nearly  energy  (and  isotope)  independent.  The 

finite  width  of  the  bands  in  SA>  BR  space  does,  however,  imply 

that  there  is  some  residual  energy  and  isotope  dependence  which 

2 

is  not  accounted  for  by  this  parameterization  of  Vqq(z) .  The 
variation  of  these  bands  over  the  range  of  8A>  6R  shown  here 
clearly  demonstrate  the  coupled  nature  of  the  dependence  of  the 
two  anisotropy  parameters  BA  and  BR»  Slight  differences  appar¬ 
ent  in  the  distribution  over  the  BA»  8R  space  between  Figures  7 
and  8  are  due  to  the  weak  variation  of  the  potential  anisotropy 
on  the  functional  form  of  VqQ(z)  as  discussed  in  section  II 
(Equation  (2.8))  . 
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Figure  9  shows  the  ratio  P  (0  -►  2)/P  (0  -*•  0)  as  a  function 

c  c 

of  incident  energy  calculated  for  H,  with  8  a  *  +0.06  [12]  and 

Z  £  L  A 

three  different  values  of  8R»  for  both  VEP  and  exp- 3  based  poten¬ 
tials.  The  experimental  ratios  PC(J  *  0  +•  2)/P(J  ■  0  -*•  0)  at  the 
four  incident  energies  of  Table  5  are  shown  as  crosses  with  as¬ 


sociated  error  bars.  Figure  10  is  the  corresponding  plot  of 
Pc(0  -+  2)/Pc(0  -+  0)  for  D 2 .  For  comparison  we  also  include  the 
results  obtained  for  a  single  8  parameterization,  as  used  in  our 
early  work  [9]  with  6  *  &A  *  6R  ■  +0 . 10 .  These  curves  are  very 

close  to  the  curves  obtained  with  8  *  -0.10,  as  expected  for 

weakly  inelastic  scattering  first  order  in  the  strength  of  the 
anisotropic  potential  component.  The  ^  results  typically  show 
threshold  resonances.  It  is  apparent  from  each  of  these  figures 
that  no  single  combination  8A»  0R  manages  to  fit  the  rotational 
transition  probability  ratios  exactly,  at  all  energies  and  for 
both  isotopes.  HoweveT,  this  empirical  fitting  procedure  does 
manage  to  distinguish  between  very  different  regions  of  the  para¬ 
meter  space  8a>  8r.  This  is  shown  in  Figure  11  where  the  ratio 

P  (0  •+  2)/P  (0  ■+  0)  is  plotted  as  a  function  of  E  for  both  VEP 

c  c  z  ^ 

and  exp-3  based  potentials  with  8*  *  +0.05,  8  °  +0.20  and  com- 

A  R 

pared  with  the  (.i)  VEP  based  potential,  8A  *  +0.06,  8R  *  +0.09 
and  (ii)  the  exp- 3  based  potential,  8A  *  +0.06,  BR  *  +0.08.  The 
values  8a  *  +0.05,  0R  *  +0.20  are  those  derived  by  Schinke  et  al . 
on  the  basis  of  the  magnitude  of  magnetic  sublevel  splittings  for 
J  *  1  H2  DSA  resonances  on  the  Ag(H0)  surface  [7,  19]  assuming  a 
VEP  based  parameterization  of  the  same  form  as  Equation  (2.8) 
with  nearly  identical  isotropic  potential  parameters  [7]. 


=  +0.05, 


It  is  clear  from  Figure  11  that  the  anisotropy  0A 
8R  *  +0.20  gives  unacceptable  values  for  the  homonuclear  rota¬ 
tional  transition  probability  ratio.  Furthermore,  whereas  for 
0A  =  +0.06,  8 ^  0.08  -  0  .09  the  ratio  Pc(0  ■+■  2)/Pc(0  •+  0)  is 
very  similar  for  both  exp- 3  and  VEP  based  potentials,  when  the 
parameters  0A  and  0R  are  significantly  different,  such  as  e.g. 

0A  =  +0.05,  0R  =  +0.20,  there  is  a  marked  difference  between  the 
energy  dependence  of  the  rotational  transition  probability  ratio 


for  the  two  different  pairs  of  functional  forms  of  vQ0  R(z)  and 


Vqq  a(z) •  Thus  one  has  to  be  very  careful  in  substantiating  the 


anisotropy  parameters  0A,  0R  obtained  from  a  fit  to  a  given 
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experimental  quantity.  When  0A  ^  0R  the  form  of  VqQ(z)  is  insen¬ 
sitive  to  the  functional  form  of  two  numerically  similar  poten¬ 


tials  Vq0(z)  (Figure  21  of  ref.  [10]).  However,  when  f*A  and  8R 


differ  by  a  large  factor,  the  form  of  VqQ(z)  becomes  sensitive 


to  the  functional  forms  of  vQq  r(z)  and  Vqq  a(z)  an<*  ^ence» 


general  sense  to  VqQ(z) . 


When  the  three  parameter  form  of  VqQCz)  was  used,  e.g.  (2.12), 
a  slight  improvement  over  the  results  of  Figures  9  and  10  was 
found.  Figure  12  shows  the  rotational  transition  probability 
ratios  for  H2  calculated  with  the  exp-3  based  potential  with 


0A  =»  +0.06,  8r  *  +0.09  and  *  0.0,  +0.05,  +0.10.  The  three 


parameter  curves  show  more  pronounced  curvature  and  thus  qualita¬ 
tively  show  increased  resemblance  to  the  experimental  results. 
Note,  however,  that  no  single  parameterization  was  found  that 
could  fit  the  experimental  ratios  for  both  homonuclear  isotopes 
at  all  collision  energies. 
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To  summarize  the  results  of  this  section,  we  find  that  a 
comparison  between  experimental  and  calculated  rotational  tran¬ 
sition  probability  ratios  for  the  weakly  inelastic  H2  and 

scattering  over  a  range  of  incident  energies  E  imply  a  range 

0  * 

of  possible  combinations  BA,  6R  for  a  given  VqQ(z)  .  In  order  to 
isolate  a  single  pair  of  values,  it  is  necessary  to  consider  addi¬ 
tional,  independent  experimental  information.  Both  the  VEP  and 
exp- 3  based  forms  of  Vpg(z,  6)  with  the  parameterization  of 
Equation  (2.8)  give  similar  bands  of  acceptable  BA>  BR  combinations, 
each  of  which  accounts  for  all  the  experimentally  measured  tran¬ 
sition  probability  ratios  to  a  first  approximation.  However,  the 
finite  width  of  these  bands  implies  some  residual  energy  dependence 
which  is  not  accounted  for  by  either  potential  with  a  two  para¬ 
meter  anisotropic  term.  This  energy  dependence  is  partially  re- 
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duced  by  the  modification  of  the  repulsive  term  in  vQ0(z)  to  in¬ 
clude  an  additional  softness  parameter. 

V.  BOUND  LEVEL  RESONANCES  AND  LINESHAPE  ANALYSIS  FOR  n-H2  AND  p-H2 

In  this  section  the  DSA  resonance  energies  for  p-H2  and  n-H2 
are  analyzed  in  order  to  derive  information  on  the  spatial  aniso¬ 
tropy  of  V00(z,  8)  which  is  complementary  to  that  obtained  via 
rotational  inelastic  scattering.  This  analysis  is  based  on  the 
fact  that  the  DSA  bound  level  resonances  for  n-H2  actually  consist 
of  multiplets,  with  this  substructure  arising  from  the  J  and  m 
dependence  of  the  resonance  energies,  En  =  en  +  en*jm‘  This 
point  was  fully  discussed  in  Section  II. B,  with  specific  reference 
being  made  to  Eqs .  (2.22-24). 
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Unfortunately,  distinct  splittings  have  not  been  observed 
in  our  DSA  data  for  either  n-H2  or  n-D2  on  Ag(lll)  [6,  10]. 

This  can  be  attributed  to  the  finite  energy  resolution  and  sen¬ 
sitivity  limitations  of  our  instrument,  and,  more  importantly, 
to  the  extremely  weak  diffractively  coupled  resonances  which  are 
observed  on  Ag(lll) ,  a  closest  packed  metallic  surface.  However, 
for  the  n  *  2  C q i 1  resonance  along  the  < 1 1 2>  direction  we  have 
observed  a  distinct  and  reproducible  example  of  resonance  line- 
shape  broadening  for  n-H2  relative  to  p-H2,  as  shown  in  Figure  13. 
The  deconvoluted  FWHM  energy  widths  of  these  resonances  are  0.19  ± 
0.04  meV  for  p-H2  and  0.35  ±  0.01  meV  for  n-H2.  These  values 
represent  the  averages  obtained  from  three  independent  measure¬ 
ments  for  each  resonance,  with  the  p-H2  and  n-H2  data  taken  sequen¬ 
tially  for  each  of  the  three  different  data  sets  in  order  to  opti¬ 
mize  the  precision  of  the  relative  lineshape  determinations.  The 
average  deconvoluted  resonance  linewidths  resulting  from  these 
measurements  are  0.47°  for  p-H2  and  0.86°  for  n-H2.  This  broaden¬ 
ing  arises  from  the  superposition  of  the  DSA  resonances  for  the 
|J,m>  =  |0,0>,  |1,0>,  and  |1,  ± 1>  substates  present  in  n-H2, 
weighted  respectively  as  **:%:>*. 

The  n  ■  2  data  for  n-H2  can  be  used  to  obtain  quantitative 
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information  on  Vqq(z)  by  carrying  out  an  accurate  computer  simu¬ 
lation  of  the  broadened  lineshape.  Such  a  computer  program,  which 
includes  a  full  convolution  over  the  experimental  parameters,  has 
been  constructed.  It  requires  as  input  the  DSA  resonance  energies 
for  each  |J,m>  state,  which  can  be  calculated  using  Eqs .  (2.23) 
and  (2.8).  It  also  requires  as  input  the  linewidth  for  each  com¬ 
ponent  of  the  multiplet.  These  were  obtained  by  assuming  that 


the  linewidth  of  each  |J,m>  sublevel  in  a  given  bound  state,  n, 
is  the  same  as  that  for  p-H2  in  that  state,  i.e.,  for  the  1 0 , 0  > 
resonance.  This  lineshape  program  has  been  used  to  examine  (i) 
for  what  range  of  anisotropy  parameters  and  a  single  broad 
resonance  is  produced  and  (ii)  which  anisotropy  parameters  give  a 
resultant  lineshape  of  0.35  meV  for  n-H2>  with  its  minimum  shifted 
by  ^<0.20  meV  from  that  for  1 0 , 0 > .  (Note  that  although  the  center- 
of-gravity  of  the  n-H2  multiplet  is  unchanged  relative  to  the  DSA 
resonance  for  p-H2,  the  lineshape  can  be  skewed  since  the  |1,  ±  1> 
peak  is  weighted  twice  that  of  the  1 1 ,0 >  resonance.) 

The  type  of  analysis  outlined  above  should  ideally  be 
carried  out  for  each  bound  level  of  Vgg(z) »  for  both  n-H2  and 
n-D2*  However,- as  was  mentioned  above,  detailed  experimental 
information  is  only  available  for  the  n  ■  2  level  of  H2»  We  shall 
therefore  restrict  our  attention  to  this  level.  Nevertheless,  it 
is  worth  mentioning  that  the  compilation  of  DSA  data  for  p-H2, 
n-H2,  o-D2  and  n-D2  presented  in  Refs.  6  and  10  exhibits  two  gen¬ 
eral  trends:  The  DSA  shifts  between  p-H2  and  n-H2  are  larger  than 
those  between  o-D2  and  n-D2,  and  the  shifts  of  both  isotopes  de¬ 
crease  as  the  bound  level  quantum  number,  n,  increases.  The  first 
part  reflects  the  different  weighting  of  the  |J,m>  states  in  n-H2 
and  n-D2,  and  the  linewidths  of  the  individual  adsorbed  |J,m> 
resonance  states.  The  second  feature  is  primarily  a  consequence 
of  the  decrease  in  expectation  value  <n|vQg(z)|n>  as  n  increases. 

The  experimental  lineshapes  of  the  p-H2  and  n-H2  resonance 
dips  were  analyzed  in  two  steps.  First,  the  deconvoluted  Lorentzian 
linewidth  of  the  p-H2  dip  was  obtained  by  fitting  the  experimental 
lineshape  with  a  forward  convolution  program.  This  program  initially 


folds  an  assumed  DSA  Lorentzian  lineshape  with  the  incident 
beam's  velocity  distribution  to  obtain  a  broadened  resonance 
lineshape  in  terms  of  resonance  energy.  It  then  transforms 
this  lineshape  into  a  function  of  incident  angle,  and  folds 
this  with  the  detector  resolution  function.  Forward  convolu¬ 
tions  of  this  type  were  then  carried  out  by  systematically  vary¬ 
ing  T,  the  HWHM  of  the  input  Lorentzian  lineshape,  until  the 
simulated  resonance  came  into  good  agreement  with  the  experi¬ 
mental  n  =  2  p-^  data.  This  comparison  was  facilitated  by 
first  fitting  the  experimental  resonance  lineshape  to  a  Gaussian 
function  after  subtracting  the  appropriate  background.  The  upper 
panel  of  Figure  14  demonstrates  that  the  experimental  p-H2  re¬ 
sonance  is  well- described  by  this  procedure,  i.e.,  by  a  Gaussian  line 
shape  with  a  broadly  sloping  quadratic  background.  The  lower 
panel  of  this  figure  shows  the  best  fit  simulation  of  this  DSA 
resonance  (solid  line),  which  occurs  for  T  =  0.05  meV.  Simulations 
using  T  *  0.5  meV  and  T  =  0.005  meV  produce  qualitatively  poorer 
fits  than  that  shown  in  Figure  16.  r  *  0.05  meV  is  used  as  the 
linewidth  for  each  of  the  multiplet  components  in  the  n-H2  simu¬ 
lation  discussed  below. 


In  the  second  step  of  the  analysis  simulated  n-H2  lineshapes 
were  constructed  by  superimposing  the  resonances  for  |J,m>  * 

1 0 , 0 > ,  |l,0>,  and  |l,  ±  1>  using  resonance  energies  calculated  for 
a  variety  of  anisotropy  parameters  6^,  8^*  These  multiplet  ener¬ 
gies  were  calculated  using  Eqs .  (2.23)  and  (2.8).  The  angular 
integrals  used  in  evaluating  these  shifts  were  <00 | P2 | 00>  *  0, 

< 10 1 P 2 1 10 >  =  +  0.4,  and  <1,  ±  l|P2|l  ±  1>  =  -  0.2.  The  expecta- 
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tion  values  <nlv0g(z)ln>  were  evaluated  for  both  the  exp- 3 


and  VEP  parameterizations  of  Vqq(z),  Eqs .  (2.9-11). 

Figure  15  displays  the  simulated  lmeshapes  for  five  sets  of 
Ba»  8r>  while  in  Figure  16  those  for  three  extreme  sets  of  ,  BR 
are  shown  along  with  the  experimental  lineshape.  The  exp- 3 
potential  was  used  for  all  of  the  simulations  shown  in  figures  15 
and  16.  Examination  of  these  two  figures  reveals  that  the  data  is 
most  closely  reproduced  by  simulations  employing  BA  ^  0.06,  with 
Br  falling  between  0.09  and  0.11.  The  lower  panel  of  Figure  16 
suggests  that  the  value  of  $R  proposed  by  Schinke,  0.20, is  too 
large  [19],  as  a  clearly  resolved  triplet  is  predicted,  in  con¬ 
tradiction  to  the  experimental  observation  of  a  single  (broadened) 
peak . 

The  results  of  the  full  BA>  0R  parameter  search  are  shown 
as  open  circles  in  Figure  8.  These  circles  denote  calculated 
lineshapes  having  energy  shifts  AE  in  the  range  0.20  ±  0.15  meV 
with  respect  to  p-H2»  and  linewidths  in  the  range  0.35  ±0.15  meV. 
It  is  interesting  to  see  that  these  circles  intersect  the  rota- 
tionally  inelastic  scattering  data  which  is  also  displayed  in 
Figure  7  at  8A  %  0.06--the  value  calculated  by  Harris  and  Feibel- 
man  [12],  and  8R  0.08--the  value  calculated  by  Harris  and  Liebsch 
for  H2/Cu(100)  [13],  but  smaller  than  the  6R  value  of  Nordlander 
and  Holmberg,  0.18  [24']. 

For  completeness,  we  end  this  section  by  reporting  our  cal¬ 
culated  shifts  for  the  |J,m>  =  |1,0>  level  (with  respect  to  |0,0>) 
for  all  of  the  bound  levels  of  H2/Ag(lll)  using  both  the  exponen¬ 
tial-3  and  VEP  parameteri zations  of  VqQ(z).  These  shifts,  for 
two  different  sets  of  8A  an<^  SR»  are  shown  in  Table  8.  The  shifts 
of  the  | J ,m>  *  | 1 ,  ±  1>  levels  can  be  easily  derived  from  Table  8 


)/2.  The  results  shown  in 


by  remembering  that  Ae  ,+1  =  -(Ae 
this  table  serve  to  emphasize  that  the  level  shifts  calculated 
with  Eq.  (2.23)  depend  not  only  on  the  values  of  SA  and  BR  employed 
in  Eq.  (2.8),  but  also  on  the  parametric  form  of  VqQ(z),  the 
laterally  averaged  isotropic  potential. 

VI .  HD  ROTATIONAL  TRANSITION  PROBABILITIES  AND  RESONANCES 
A.  Experimental  Results 

In  this  section  the  experimentally  determined  rotationally  in¬ 
elastic  transition  probabilities  for  HD  are  presented.  These  an¬ 
gular  distributions  were  taken  at  incident  angle  increments  of  5°, 
spanning  the  range  from  0^  =  25°  to  0.  =  85°,  for  two  beam  energies 
Figure  17  shows  a  typical  angular  distribution  and  the  correspond¬ 
ing  Ewald  representation  of  the  scattering  kinematics.  All  of  the 
diffractive  and  rotationally  inelastic  peaks,  labelled  as  J ^ 
in  this  figure,  are  directly  interpretable  in  terms  of  (i)  paral¬ 
lel  momentum  conservation,  Eq.  (3.1),  and  (ii)  conservation  of 
total  energy,  Eq.  (3.2).  A  full  discussion  of  the  scattering  kine¬ 
matics  for  HD/Ag(lll)  has  been  previously  presented  [9],  and  will 
not  be  repeated  here.  The  main  features  of  the  HD/Ag(lll)  system 
are  the  extremely  high  rotationally  inelastic  transition  probabil¬ 
ities  for  HD  relative  to  H2  and  D2  (see  Figs.  2,  3,  and  17),  and 
the  presence  of  only  weak  diffraction  for  G  f  0.  The  first  obser- 
vation  can  be  attributed  to  the  offset  of  the  HD  center-of-mass 
from  its  geometric  center,  while  the  second  indicates  that  the 
HD/Ag(lll)  potential  has  very  low  corrugation.  This  latter  point 
was  discussed  in  Section  IV  and  allows  us  to  treat  Ag(lll)  to 
first  approximation  as  a  flat  surface  in  the  scattering  calcula¬ 
tions  presented  in  this  paper.  Similar  observations  have  also 
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been  recently  reported  for  HD/Pt(lll)  [22J.  It  is  apparent 
from  Figs.  4  and  17  that  J  s  1  or  higher  rotational  states  are 
nearly  absent  in  the  incident  beam,  as  evidenced  by  the  small 
AJ  *  1  -*■  2  peaks  in  these  figures.  We  shall  therefore  treat 
the  incident  HD  |J,m>  ensemble  as  consisting  exclusively  of  the 
1 0 , 0 >  state  in  our  data  analysis  and  theoretical  calculations. 
Finally,  our  rotational  assignments  (made  initially  on  the  basis 
scattering  kinematics)  have  been  independently  confirmed  by 
recording  the  actual  time-of-flight  distributions  of  the  scat¬ 
tered  molecules,  as  shown  in  Figure  18  for  0^  “  25°.  These  four 
velocity  distributions  are  in  excellent  agreement  with  those 
predicted  by  the  conservation  of  parallel  momentum  and  total 
energy. 

Six  HD  angular  distributions ,  taken  at  incident  angle  in¬ 
crements  of  10°,  are  shown  in  Figure  19.  These  distributions 
demonstrate  that  additional  inelastic  channels  become  accessible 
as  6^  decreases  towards  the  surface- normal,  i.e.,  as  the  normal 
component  of  the  incident  beam  energy  increases.  The  deconvolu- 
ted  inelastic  transition  probabilities  for  the  room  temperature 
and  cryogenically  cooled  HD  beams  are  given  in  Tables  9  and  10, 
respectively.  These  results  were  derived  from  the  experimental 
angular  distributions  using  the  same  procedures  as  were  used  in 
Section  IV  for  treating  H2  and  D2.  The  theoretical  calculations 
presented  in  the  next  section  were  compared  to  the  Debye-Waller 
corrected  data  listed  in  Table  9  after  rescaling  the  transition 
probabilities  such  that  they  sum  to  unity  for  each  angular  dis¬ 
tribution.  This  normalization  procedure  ignores  the  peaks  asso¬ 
ciated  with  G  f  0 ,  which  account  for  less  than  51  of  the  reflected 


flux.  Our  data  analysis  and  normalization  procedures  were  checked 

by  comparing  the  transition  probabilities  derived  from  Tables  9 

and  10  as  a  function  of  E  .  These  energy  dependent  final  state 

zi 

distributions  (for  two  different  beam  energies)  were  found  to  be 
superimposable- -validating  our  analysis  procedures,  and  again  con¬ 
firming  that  the  flat  surface  approximation  is  valid  for  the 
HD/Ag(lll)  system. 

B.  Theoretical  Analysis 

In  this  section  we  compare  the  theoretical  predictions  for 
the  HD  rotational  transition  probabilities  and  RMSA  resonances 
with  the  experimental  results,  using  the  anisotropy  parameters 
that  were  derived  in  sections  IV  and  V  from  the  H 2  and  D2  results. 

As  discussed  earlier  in  section  II,  both  the  rotationally  inelas¬ 
tic  scattering  and  resonance  energies  for  HD  are  relatively  in¬ 
sensitive  to  small  changes  in  the  anisotropy  of  the  potential, 

Vqq(z),  and  we  include  these  results  here  primarily  for  confirma¬ 
tion  of  the  conclusions  from  and  D2,  to  differentiate  between 
widely  separated  areas  of  the  BR  parameter  space, and  for  com¬ 
pleteness.  We  note  that  since  the  isotropic  potential  was  deter¬ 
mined  by  the  DSA  resonances  of  H2  and  D2  for  J  =  0,  its  determina¬ 
tion  is  independent  of  the  values  of  B^  arid  BR  assumed.  Therefore 
no  iteration  of  the  isotropic  and  anisotropic  potentials  is  required. 
The  calculated  rotational  transition  probabilities  Pc(0  -+  J'), 

J’  =  0,1, 2, 3  are  shown  in  Figures  20  and  21  as  a  function  of  E  , 
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for  the  VEP  based  potential  with  =  0.06,  6R  =  +0.09  and  for  the 
exponential- 3  based  potential  with  $A  =  +0.06,  8R  *  +0.08,  respect¬ 
ively.  The  experimental  values  (unitarized  transition  probabilities 


in  these  figures.  The  results 
similar,  and  are  in  good  agree- 
.  Although  the  HD 


transition  probabilities  show  little  sensitivity  to  email  changes 
in  the  anisotropy  parameters  6A  and  BR,  they  are  qualitatively 
affected  by  relatively  large  changes  in  BA  and  gR.  For  example, 
in  Figure  22  we  show  the  calculated  transition  probabilities  for 
the  exponential- 3  potential  with  8A  *  +0.05,  8R  =  +0.20.  This 
gives  significantly  poorer  agreement  with  the  experimental  prob¬ 
abilities  than  8a  *  +0.06,  BR  «  +0.08  (Tig-  21).  When  the  VEP 
parameterization  is  used  with  $A  =  +0.05,  BR  =  +0.20  the  overall 
quality  of  the  fit  is  also  poorer  than  the  results  obtained  with 
8a  =  +0.06,  6r  a  +0.09  (Fig.  20),  but  not  to  the  striking  extent 
that  was  seen  between  Figs.  21  and  22.  We  therefore  see  that 
even  for  HD,  where  the  rotational  inelasticity  is  largely  dominated 
by  the  high  mass  asymmetry  of  the  molecule,  the  energy  dependence 
of  the  inelasticity  can  clearly  distinguish  between  acceptable 
and  unacceptable  regions  of  the  $A»  8R  parameter  space. 

Comparison  of  the  calculated  RMSA  resonance  energies  with  the 
experimentally  measured  values  [9]  is  made  in  Table  11.  The  experi¬ 
mental  values  derived  from  the  VEP  potential,  Ee,  are  compared 
here  to  the  zeroth  order  BVRR  energies  derived  from  the  VEP  poten¬ 
tial,  E^,  and  to  the  exact  resonance  energies  Ec*  derived  from 
six  parameterizations  i  of  the  potential  VQ0(z,  e) :  (1)  VEP; 

8a  *  +0.06,  0R  •  +0.09  (2)  exp-3;  8A  -  +0.06,  BR  *  +0.08  (3)  exp-3; 

6A  -  +0.06,  gR  -  +0.09,  a2  •  +0.10  (4)  VEP;  $A  -  +0.05,  0R  -  +0.20 

(5)  VEP;  BA  *  8R  *  +0.10  (6)  VEP;  BA  *  6R  ■  -0.10.  The  results 

from  potentials  (1),  (2),  and  (3)  which  all  have  anisotropic  para¬ 
meterization  in  the  region  of  our  previous  conclusions  derived  from 
H2  and  D2,  differ  by  at  most  0.2  meV,  which  is  the  magnitude  of  the 


experimental  precision.  Comparing  to  the  experimental  resonance 
energies,  with  the  exception  of  the  three  highest  energy  resonances 
for  adsorbed  J  =  3,  the  calculated  energies  fall  within  0.2  meV  or 
less  of  the  experimental  values.  The  J  =  3,  n  =  3,4  resonances 
are  predicted  to  be  too  low  by  several  meV  for  all  parameterizations 
shown;  this  was  seen  for  all  combinations  g^,  g^  probed  in  figures 
7  and  8.  Comparing  now  with  results  from  potential  (4)  (BA  *  +0.05, 
0R  =  +0.20),  we  see  very  similar  behavior.  As  discussed  in  Section 
II,  the  RMSA  resonance  energies  are  strongly  dominated  by  the  ani¬ 
sotropy  arising  from  the  molecular  asymmetry.  In  contrast  to  the 
energy  dependence  of  the  HD  rotational  transition  probabilities, 
it  is  not  possible  to  distinguish  between  these  two  regions  of 
$A»  0R  parameter  space  from  the  RMSA  resonance  positions. 

As  a  final  comparison  we  have  also  included  in  Table  11  the 
resonance  energies  predicted  by  the  single  g  parameterization 
(3  =  $A  =  Br)  *  potentials  (5)  and  (6).  These  show  that  both  the 
sign  and  magnitude  of  the  shift  depend  on  the  sign  of  g .  It  is 
interesting  that  the  value  g  =  -0.10  gives  approximately  the  same 
shifts  as  8a  =  +0.05,  g^  ■  +0.20,  an  observation  which  can  be 
rationalized  in  terms  of  the  relative  sizes  of  the  expectation 
values  of  vQq  r(z)  and  Vqq  a(z)  [25].  This  equivalence  of  the 
RMSA  shifts  is  to  be  contrasted  with  the  very  different  behavior 
seen  for  the  ^  rotational  transition  probability  ratios  for  these 
two  anisotropic  parameterizations,  within  the  VEP  based  form  of 
Yq0(z,  6).  Such  an  observation  demonstrates  once  again  the  need 
to  analyze  both  transition  probabilities  and  resonance  energies 
in  order  to  extract  unambiguous  information  about  the  anisotropic 
parameters . 


VII.  SUMMARY  AND  CONCLUSION 


In  this  paper  we  have  investigated  the  anisotropic  component 
of  the  laterally  averaged  interaction  potential,  VQ0(z,  e) ,  for 
molecular  hydrogen/Ag(lll) ,  by  employing  a  combination  of  experi¬ 
mental  measurements  and  theoretical  calculations.  The  calculations 
utilized  a  laterally  averaged  potential  which  was  expanded  in  a 
Legendre  series  for  a  symmetric  (homonuclear)  system,  Eqs .  (2.5-6), 

in  which  the  anisotropic  component  was  given  by  a  two  parameter 
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model,  Eq.  (2.8).  This  two  parameter  representation  of  vQ0(z) 

allowed  us  to  systematically  examine  the  orientation  dependence 

of  both  the  attractive  and  repulsive  parts  of  the  potential.  A 
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three  parameter  model  of  Vqq(z)  was  also  briefly  examined  which 

2 

allowed  the  softness  of  the  repulsive  component  of  vQ0(z)  to  be 
independently  varied.  The  isotropic  potential  used  in  these  cal¬ 
culations  was  represented  by  variable  exponent  and  exponential- 3 
parameterizations ,  Eqs.  (2.9-11),  that  were  previously  shown  to 
reproduce  the  experimental  bound  state  resonance  spectra  for  p-I^ 
and  0-D2  on  Ag(lll)  [6,10].  The  exponential- 3  potential  was  used 
to  a  greater  extent  in  our  analysis  than  the  variable  exponent 
potential  as  it  has  the  proper  form  predicted  by  ab  initio  calcu¬ 
lations  . 

Our  analysis  began  by  comparing  experimentally  derived  rota- 
tionally  inelastic  transition  probabilities,  taken  as  a  function 
of  collision  energy  normal  to  the  surface,  with  those  resulting 
from  close  coupled  quantum  scattering  calculations.  Extensive 
iterations  of  the  anisotropy  parameters  8A  and  8R  in  these  calcu¬ 
lations  revealed  that  the  energy  dependent  H2  and  D2  transition 


probabilities  fall  along  a  common  band  of  BA,  BR  values  (Figs.  7 
and  8),  with  small  BR  correlating  with  small  SA-  Comparable 
results  were  found  using  the  VEP  and  exponential- 3  isotropic 
potentials,  with  the  exponential-3  results  showing  slightly 
better  agreement  with  the  transition  probabilities  as  a  function 
of  energy.  When  0A  is  set  at  the  theoretical  value  predicted  by 
Harris  and  Feibelman  [12],  0A  *  0.06,  a  0R  value  of  0.09  is  found 
for  the  VEP  parameterization  of  vjjgCz),  an<*  for  the 

exponential-3  based  potential.  These  values  for  0R  are  in  good 
agreement  with  the  ab  initio  result  for  H2/Cu(100)  given  by  Harris 
and  Liebsch  [13],  BR  ^  0.08,  but  are  smaller  than  the  value  calcu¬ 
lated  by  Norlander  and  Holmberg,  BR  'v  0.18  [24].  Inclusion  of  an 
additional  exponentially  decaying  term  in  the  repulsive  anisotropic 
potential  causes  some  improvement  between  the  energy  dependence  of 
the  experimental  and  calculated  probabilities  to  be  realized.  For 
example,  an  exponential- 3  based  potential  which  includes  this  addi¬ 
tional  softness  parameter  gives  a  better  fit  to  the  data  with 
Ba  *  0.06,  0R  ■  0.09,  and  a2  *  0;10  than  was  originally  obtained 
with  Ba  *  0.06,  0R  -  0.08,  and  a2  *  0.00.  We  also  note  that  the 
transition  probabilities  for  the  homonuclear  hydrogen  isotopes 
are  quite  sensitive  to  changes  in  0A  and  0R,  and  can  therefore  be 
used  to  distinguish  between  different  regions  of  the  0A>  8R  para¬ 
meter  space.  This  was  demonstrated  in  Fig.  11,  where  the  H2  rota¬ 
tional  transition  probabilities  calculated  with  0R  'v  0.20  were 
found  to  be  in  poorer  agreement  with  the  data  than  those  calculated 
using  0R  'v  0.09  (especially  for  the  exponential- 3  potential). 
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The  DSA  resonance  energies  for  p-H2  and  n-H2  were  then 
analyzed  in  order  to  examine  how  the  magnetic  sublevel  splittings 
for  J  ■  1  n-H2  evolve  as  a  function  of  BA»  BR  and  bound  level  n. 
Unfortunately,  distinct  splittings  have  not  been  observed  in  our 
DSA  data  for  either  n-H2  or  n-D2  on  Ag(lll)  [6,10].  This  can  be 
attributed  to  the  extremely  weak  nature  of  the  diffractively  coupled 
resonances  which  are  observed  on  Ag(lll) ,  a  closest  packed  metallic 
surface,  and  to  the  finite  energy  resolution  of  our  instrument. 
However,  the  n  *  2  (g®)  resonance  for  n-H2  along  the  <112>  direc¬ 
tion  exhibited  reproducible  linewidth  broadening  relative  to  p-H2, 
and  was  examined  using  a  detailed  forward  convolution  procedure. 

The  results  of  this  analysis,  given  as  a  function  of  BA  and  BR  in 
Fig.  8,  intersect  the  rotationally  inelastic  H2  data  at  8A  ^  0.06, 

8r  ^  0.08.  Although  the  quantitative  information  derived  from  this 
lineshape  analysis;  pertaining  to  the  anisotropic  potential  is  the 
least  definitive  evidence  presented  in  this  paper  (due  to  extensive 

convolutions  and  low  signal  level),  it  does  serve  to  demonstrate 

2 

that  this  type  of  experiment  can  be  used  to  examine  Vqq(z)  in  a 
manner  complementary  to  inelastic  scattering.  However,  we  wish 
to  strongly  emphasize  that  the  calculated  energy  level  shifts  for 
the  |J,m>  -  |1,0>  and|l,  ±  1>  states  relative  to  |0,0>  depend  not 

only  on  the  values  of  BA  and  BR»  i.e.,  on  the  anisotropic  potential, 
but  also  on  the  parametric  form  of  Vqq(z),  the  laterally  averaged 
isotropic  potential.  This  is  demonstrated  in  Table  8,  where  the 
shifts  for  J  =  1  n-H2  exhibit  an  obvious  dependence  on  the  isotropic 
potential  used  in  the  calculation. 

Finally,  the  BA  and  8R  values  determined  in  the  H2  and  D2 
experiments  were  used  to  calculate  rotational  transition  probabilities 
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and  RMSA  resonance  positions  for  HD  incident  upon  Ag(lll)  modelled 
as  a  flat  surface.  These  calculated  values  were  then  compared  to 
our  experimental  results.  We  find  that  although  the  rotationally 
inelastic  transition  probabilities  and  RMSA  resonance  energies 
for  HD  are  relatively  insensitive  to  small  changes  in  the  aniso¬ 
tropy  of  the  potential,  the  results  tend  to  support  the  conclusions 
drawn  from  the  H 2  and  D2  experiments.  In  particular,  the  rotation- 
ally  inelastic  scattering  results  for  HD  differentiate  between  the 
Ba  *  0.06,  Br  =  0.09  and  BA  =  0.05,  BR  =  0.20  parameter  sets,  with 
the  former  giving  substantially  better  agreement  to  the  data. 

This  can  be  seen  for  the  exponential- 3  based  potential  in  Figs.  21 
and  22.  In  contrast  to  this,  the  HD  RMSA  resonance  energies  are 
found  to  exhibit  very  similar  J- dependent  level  shifts  using  the 
Ba  *  0.06,  Br  *  0.09  and  BA  *  0.05,  8^  *  0.20  parameter  sets.  This 
surprising  insensitivity  to  the  anisotropic  potential  arises  for 
the  heteronuclear  isotope  since  the  RMSA  resonance  energies  are 
strongly  dominated  by  the  molecular  asymmetry  of  this  system. 

In  conclusion,  we  believe  this  investigation  has  demonstrated 
that  in  order  to  make  an  accurate  and  unambiguous  determination  of 
a  molecule-surface  potential,  a  wide  combination  of  complementary 
measurements  and  calculations  are  needed.  Diffractive  selective 
adsorption  and  thermal  desorption  measurements  can  be  used  to 
examine  the  orientation  independent  part  of  the  potential,  as  was 
done  in  our  first  paper  of  this  series  on  hydrogen/Ag(lll)  [10]. 
Diffraction  and  diffraction-rotation  measurements  can  probe  the 
periodic  part  of  the  potential,  as  will  be  discussed  in  a  future 
publication  [11].  Finally,  as  demonstrated  in  this  paper,  rota¬ 
tionally  inelastic  scattering,  and  selective  adsorption  resonances 
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which  exhibit  J  and  m  dependent  level  shifts,  can  be  used  to  examine 
the  anisotropic  part  of  the  potential.  Using  this  combination  of 
experimental  measurements  and  calculations  we  have  determined 
attractive  and  repulsive  anisotropy  parameters  which  agree  well 
with  calculated  values  [12,13].  This  study,  together  with  the 
previous  one  dealing  with  the  isotropic  potential,  provides  the 
most  complete  low  energy  molecule- surface  potential  determined  to 
date . 
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TABLE  1 

EMPIRICAL  ROTOR  ENERGIES,  E(J)  ,  ^  FOR  HD  USED  IN  ALL 
CLOSE  COUPLED  SCATTERING  CALCULATIONS.  RIGID  ROTOR 
ENERGIES  ARE  SHOWN  FOR  COMPARISON.  THE  GROUND  VIBRATIONAL 
STATE  OF  HD  IS  ASSUMED.  ALL  ENERGIES  ARE  IN  meV. 


Rotor  State  J 

Empirical  E(J) 

Rigid  Rotor  Energy 

0 

0.000 

0.000 

1 

11.063 

11.320 

2 

33.114 

33.961 

3 

66.000 

67.922 

4 

109.502 

113.204 

5 

163.331 

169.806 

(a)  H.  W.  Woolley, 

R.  B.  Scott,  and  F. 

G.  Brickwedde, 

J.  Res.  Nat.  Bur.  Stand.  44,  379  (1948). 
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TABLE  2 

PARAMETERS  FOR  THE  ISOTROPIC  COMPONENT  v°Q(z)  OF  THE 
LATERALLY  AVERAGED  INTERACTION  POTENTIAL 


Potential  Parameters 

Variable  exponent  D  =  31.535  meV 

X  =  0.6335  a^ 

P  =  4.2925 

2  =  2.3441  an 

e  0 

Exponential- 3  D  =  32.4616  meV 

a  =  1.2789  aQ~ 

zg  =  3.7572  a^ 


TABLE  5 


ROTATIONALLY  INELASTIC  TRANSITION  PROBABILITIES  FOR 
n-H2  AND  p-H2.  Tg  =  125  K,  Ei  =  77.5  meV  (n-H2) ,  Et  *  79.6  meV 
(p-H2) .  NOTE  THAT  THE  EXPERIMENTAL  SPECULAR  TRANSITION 
PROBABILITIES  INCLUDE  SEVERAL  J  *  n  +  n  COMPONENTS. 


““ 

Gas 

0. 

i 

Spec. 

0  +  2 

o 

+ 

o 

0  +  2 

(0  + 

2)/(0  +  0) 

25 

.224 

.00182 

.207 

.00598 

.0289 

±  .0043 

30 

.207 

.00151 

.180 

.00456 

.0254 

±  .0038 

n_H2 

35 

.214 

.00082 

.172 

.00228 

.0133 

±  .0027 

40 

.207 

.00044 

.154 

.00111 

.0072 

±  .0014 

p-H2  27.5  .210 


00510 


448 


01025 


0229  ±  .0029 
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TABLE  6 

ROTATIONALLY  INELASTIC  TRANSITION  PROBABILITIES  FOR 
n-D2.  Tg  *  125  K,  Ei  -  82.6  meV.  NOTE  THAT  THE 
EXPERIMENTAL  SPECULAR  TRANSITION  PROBABILITIES  INCLUDE 
SEVERAL  J  =  n  ♦  n  COMPONENTS. 


Experimental 


Debye-Waller  Corrected 


0. 

X 

Spec . 

0  -  2 

0  +  0 

0  +  2 

(0  +  2)/(0  +  0) 

25 

.0340 

.00195 

.202 

.0203 

.1005  ±  .0231 

.0447 

.00214 

.234 

.0194 

.0829  ±  .0149 

35 

.0547 

.246 

.0164 

.0668  ft  .0120 

.0650 

.00184 

.253 

.0123 

.0486  ±  .0063 

45 

.0868 

.00187 

.288 

.0106 

.0369  ±  .0048 

.1102 

.00178 

.311 

.00B5 

.0274  ft  .0036 
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®The  perpendicular  energies  Ez  ■  cos^0j  correspond  to  the  experimental  energies.  These 
values  of  j}  u  were  obtained  by  comparison  of  the  calculated  ratio  of  transition  probabilities  P(J-0-*2)/ 
P(J-0->0)  to  the  experimental  ratios.  Both  £.  and  ^  are  incremented  by  0.01  In  this  search.  Inter-1 
polated  values  are  given  where  appropriate  ana  are  bracketed  to  indicate  this. 
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TABLE  8 

CALCULATED  ENERGY  LEVEL  SHIFTS  Aen  FOR 
J  -  1,  m  -  0  HYDROGEN  ADSORBED  ON  Ag(lll) 
FOR  THE  EXPONENTIAL- 3  AND  VARIABLE  EXPONENT 
PARAMETERIZATIONS  OF  VQ0(z)(a) 
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TABLE  9 


ROTATIONALLY  INELASTIC  TRANSITION  PROBABILITIES  FOR 
HD/Ag (111) .  T  =  155  K,  E.  =  109.7  meV 

3  X 


Experimental 


ei 

0  -  0 

0-1 

0-2 

0-3 

0-0 

0-1 

0-2 

0  - 

25 

.0115 

.0326 

.0455 

.00665 

.217 

.539 

.571 

.054 

30 

.0097 

.0412 

.0471 

.00440 

.152 

.563 

.483 

.054 

35 

.0104 

.0487 

.0512 

.00289 

.131 

.539 

.430 

.015 

40 

.0104 

.0649 

.0593 

— 

.105 

.574 

.397 

— 

45 

.0118 

.0925 

.0614 

— 

.095 

.649 

.325 

— 

50 

.0142 

.1217 

.0651 

— 

.090 

.678 

.273 

■  -  -  - 

55 

.0212 

.1638 

.0509 

— 

.108 

.729 

.170 

— 

60 

.0338 

.2186 

— 

— 

.139 

.789 

— 

— 

65 

.0546 

.2235 

— 

— 

.186 

.666 

— 

— 

70 

.0737 

.2121 

— 

— 

.214 

.537 

— 

— 

75 

.2472 

— 

— 

— 

.691* 

— 

— 

— 

80 

.2405 

— 

— 

— 

.774* 

— 

— 

— 

85 

.2062 

— 

— 

— 

.846* 

— 

— 

— 

Debye-Waller  Corrected 


♦Corrected  for  the  geometrical  loss  of  incident  beam  flux  which 


occurs  for  9^  >  75  . 


f. (.'f.  -  . -  .  -  .  <■  ••  V.  .■  .'  V  -*  ^  ’.'j  *  .  .  '  V.S  .s 


TABLE  10 

ROTATIONALLY  INELASTIC  TRANSITION  PROBABILITIES  FOR 
HD/Ag(lll) .  Tg  *  155  K,  E.^  =  47.6  meV 


- - 

Experimental 

Debye-Waller  Corrected 

e. 

i 

0-^0 

0  -*■  1 

0  -*■  2 

0  -*•  0 

0  -►  1  0  -►  2 

29.2 

.0066 

.0546 

.0147 

.0337 

.244  .0492 

40.5 

.2830. 

.0125 

.0788 

.0516 

.285 

58.6 

.0352 

;  .1009 

— 

.1020 

.255 

74.7 


.1186 


2730 


COMPARISON  OF  RMSA  RESONANCE  ENERGIES  FOR  HD  OBTAINED  BY  CALCULATION  E  WITH  THE 
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FIGURE  CAPTIONS 

Figure  1  Coordinate  system  for  a  diatomic  molecule  AB  incident 


Figure  2 


Figure  3 


Figure  4 


upon  a  stationary  corrugated  surface.  0  is  the  origin 
of  coordinates,  fixed  at  some  reference  point  in  the 
surface  unit  cell,  C  the  midpoint  of  the  internuclear 
distance  AB  (i.e.,  the  center  of  geometry),  E  the 
molecular  center  of  mass.  PQ  is  the  projection  of 
AB  in  the  surface  xy  plane,  S  and  T  are  the  projec¬ 
tions  of  C  and  E,  respectively.  The  locations  of  C 
and  E  in  the  xyz  frame  are  thus  given  by  vectors 
(R* ,z')  and  (R,z),  respectively,  where  R  *  (x,y) . 

Atoms  A  and  B  have  masses  m^  and  m2  and  coordinates 
(R2,z1)  and  (R2,z2),  respectively.  The  internal 
degrees  of  freedom  are  expressed  in  the  spherical 
polar  coordinates  (r,8,4>)  where  r  is  the  internuclear 
distance  AB,  and  0  and  <)»  the  polar  and  azimuthal 
angles,  measured  in  the  xyz  frame.  EC,  the  offset 
of  the  center  of  mass  from  the  center  of  geometry, 
is  denoted  by  S  and  is  equal  in  magnitude  to  rCm^-rn^/ 
2(m1+m2) . 

Four  angular  distributions  for  n-H2/Ag(lll)  with 
4>  *  15°  and  Tg  *  125  K. 

Four  angular  distributions  for  n-D2/Ag(lll)  with 

4>  =  15°  and  T  *  125  K. 
s 

Comparison  of  HD  angular  distributions  taken  at  dif¬ 
ferent  beam  energies  but  with  the  same  E  . 

Z  4 


Figure 


Computer  fit  used  to  deconvolute  the  J  =  0  -*■  2  in¬ 


elastic  transition  probability  for  D2/Ag(lll)  at 
8i  =  45°.  The  solid  line  is  the  fit  to  the  data 
(X)  while  the  broken  line  is  the  fit  to  the  diffuse 
background.  The  computer  model  clearly  reproduces 
the  experimental  data  very  well. 

Figure  6  Rotational  transition  probability  ratio  P(J  =  0  -*■  2)/ 
P(J  =  0  -*■  0)  for  H2  measured  as  a  function  of  BR  for 
a  range  of  values  of  BA>  at  incident  energy  Ez  =  53.2 
meV.  Both  parameters  6A  and  BR  are  incremented  by 
0.01.  The  experimentally  measured  ratio  is  given  by 
the  solid  horizontal  line  and  the  associated  error 
bounds  shown  as  dashed  lines. 

Figure  7  Anisotropic  parameter  space  BA,  BR  for  the  VEP  poten¬ 
tial  showing  the  interaction  between  regions  by  fit¬ 
ting  H2  and  D2  rotationally  inelastic  transition 
probabilities.  Crosses (x)  denote  combinations  de¬ 
rived  from  the  H2  rotationally  inelastic  transition 
probabilities,  while  pluses (+)  are  for  D2.  The  band 
of  B^  values  for  given  BA  reflects  the  energy  depend¬ 
ence  apparent  in  Table  7. 
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Figure  8  Anisotropic  parameter  space  8A,  $R  for  the  exp- 3 
potential  showing  the  interaction  between  regions 
derived  by  fitting  (i)  H2  and  D2  rotationally  in¬ 
elastic  transition  probabilities,  and  (ii)  a  single 
|Jm>  -  averaged  first  order  shift  AEn  and  linewidth 
broadening  for  n-H2  adsorbed  in  the  n  *  2  level. 
Crosses (x)  denote  combinations  derived  from  the  H2 
rotationally  inelastic  transition  probabilities 
while  pluses (+)  are  for  D2.  Circles(o)  denote  com¬ 
binations  derived  from  the  n-H2  shift  and  linewidth 
broadening  in  the  n  =  2  level  (see  Section  V) .  The 
linewidth  is  measured  as  the  FWHM  of  the  n-H2  distri¬ 
bution.  Only  distributions  with  one  resultant  peak 
or  two  peaks  separated  by  less  than  0.5  meV  are  in¬ 
cluded.  Further  selection  is  made  by  requiring  the 
averaged  shift  AE  to  be  in  the  range  +0.20  *  0.15  meV 
and  the  linewidth  to  be  in  the  range  0.35  ±  0.15  meV. 

Figure  9  Rotational  transition  probability  ratio  P(J  *  0  •+  2)/ 

P(J  *  0  -*•  0)  for  H2  as  a  function  of  Ez,  calculated 

0  2 

for  different  forms  of  VqQ(z)  and  vQ0(z).  Results 
from  the  VEP  based  potential  Equation  (2.9)  with 
r  *  +0.08,  +0.09,  +0.10  are  shown  as 

solid  lines  ( - ).  Results  for  the  exponential-3 

based  potential,  Equation  (2.10)  with  BA  *  +0.06, 

Br  *  +0.07,  +0.08,  +0.09  are  shown  as  dashed  lines 


'*  -'-.*VVV's 
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( —  —  — ) .  Results  for  the  VEP  based  potential 
with  the  single  $  anisotropy  B  =  $A  =  BR  =  +0.10 
are  shown  as  dot-dashed  lines  ( —  •  —  •  — ).  Cal¬ 
culations  were  performed  at  increments  of  2.0meV  in 
Ez.  Experimental  values  are  denoted  by  crosses(x) 
and  associated  error  bars.  The  data  point  denoted 
by  a  dot  (•)  was  determined  using  a  p-f^  incident 
beam. 

Figure  10  Rotational  transition  probability  ratio  P(J  =  0  •+  2)/ 
P(J  =  0  -*•  0)  for  as  a  function  of  E  ,  calculated 
for  different  forms  of  vQ0(z)  and  vQ0(z).  Results 
from  the  VEP  based  potential  Equation  (2.9)  with 
6a  *  +0.06,  R  =  +0.08,  +0.09,  +0.10  are  shown  as 

solid  lines  ( - ) .  Results  for  the  exponential- 3 

based  potential.  Equation  (2.10)  with  BA  *  +0.06, 

BR  =  +0.07,  +0.08,  +0.09  are  shown  as  dashed  lines 

( - ).  Results  for.  the  VEP  based  potential  with 

the  single  anisotropy  B  =  BA  =  8R  =  +0.10  are 
shown  as  dot-dashed  lines  ( —  •  —  •  — ) .  Calcu¬ 
lations  were  performed  at  increments  of  2.0meV  in 
Ez.  Experimental  values  are  denoted  by  crosses (x) 
and  associated  error  bars. 

Figure  11  Rotational  transition  probability  ratio  P(J  =  0  -+  2)  / 

P(J  =  0-+0)  for  H2  as  a  function  of  Ez  showing  the 

effects  of  large  variations  in  the  parameterization 
2 

of  Vqq(z).  The  solid  line  ( - )  is  the  result  for 

the  exponential-3  based  potential  with  parameters  of 


62 


of  Table  2  and  anisotropy  BA  =  +0.06,  8R  =  +0.08 
while  the  ( —  •  —  •)  line  is  the  result  for  the 
exponential-3  potential  with  anisotropy  BA  =  +0.05, 

Br  =  +0.20.  The  ( —  ••  —  ••)  line  is  the  result 
for  the  VEP  based  potential  with  parameters  of 
Table  2  and  anisotropy  BA  =  +0.06,  8R  =  +0.09  while 

the  dashed  line  ( —  - )  is  the  result  for  the  VEP 

potential  with  anisotropy  BA  =  +0.05,  BR  =  +0.20. 
Calculations  were  performed  at  increments  of  2.0 
meV.  Experiment  values  are  donated  by  crosses (x) 
for  n-H2  derived  values  and  a  dot  (*)  for  the  p-H2 
derived  value. 

Figure  12  Rotational  transition  probability  ratios  P(J  =  0  -*•  2)/ 
P(J  =  0  -+  0)  for  H2  as  a  function  of  Ez  for  the  expon¬ 
ential-3  based  potential  with  the  three  parameter 
anisotropic  term  given  by  Equation  (2.12).  Results 
are  shown  for  anisotropy  parameters  $A  =  +0.06, 

8r  =  +0.09  with  $2  =  +0.05,  +0.10  ( - ).  The  iso¬ 

tropic  parameters  are  given  in  Table  2.  Calculations 
are  carried  out  on  a  grid  of  2.0meV  in  Ez<  Experi¬ 
mental  values  are  denoted  by  crosses(x)  and  associa¬ 
ted  error  bars.  The  data  point  denoted  by  a  dot  (•) 
was  determined  using  a  p-H2  incident  beam. 


Figure  13  Selective  adsorption  linewidth  broadening  for  n-H2 
relative  to  p-H2<  The  total  shift  in  dip  angle  is 
due  to  both  small  |j,m>  -dependent  level  shifts  for 
n-H2  and  the  difference  in  beam  energies. 

Figure  14  Selective  adsorption  resonance  lineshape  for  the  n  =  2 
level  of  the  p-H2/Ag(lll)  physisorption  potential. 

The  upper  panel  shows  the  best  Gaussian  fit  (solid 
line)  to  the  data  (•)  including  the  sloping  back¬ 
ground  (see  Fig.  IS) .  The  lower  panel  contains  the 
same  Gaussian  fit  to  the  data  as  above  but  without 
background  ( —  —  — )  and  the  simulated  lineshape 

( - •)  •  The  computer  simulation  is  a  full  convolu¬ 

tion  of  a  Lorentzian  lineshape  (HWHM  0.05  meV)  with 
the  detector  resolution  function. and  the  velocity 
spread. of  the  incident  beam. 

Figure  15  n-H2  selective  adsorption  lineshape  simulations  for 

the  n  =  2  resonance  shown  in  Fig.  15  as  a  function  of 
Ba,  Br.  Input:  r  *  0.05  meV,  beam  velocity  spread 
0.09  meV,  detector  resolution  0.66°. 

Figure  16  Three  additional  simulations  of  the  n  =  2  resonance  at 
(j>  =  0°  (see  Figure  15)  for  a  wide  parameter  sampling. 
Dashed  curves:  experimental  lineshape;  solid  lines: 
computer  simulations  including  the  instrument  transfer 
function. 
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Figure  17  Angular  distribution  of  the  reflected  number 

density  for  HD  scattering  from  Ag(lll)  along  the 
<112>  direction.  Rotationally  inelastic  transitions 
are  labelled  -*■  J^.  Vertical  lines  denote  the 
predicted  peak  locations  for  a  monoenergetic  beam 
with  k^  =  12.7  A  The  lower  section  is  an  Ewald 
representation  of  the  displayed  angular  distribution. 


Figure  18  HD  time-of- flight  distributions  taken  at  the  peak  posi¬ 
tions  of  the  various  elastic  and  rotationally  inelastic 
channels . 


Figure  19  Angular  distributions  of  the  reflected  number  density 

for  HD/Ag(lll)  taken  at  six  incident  angles.  Rotation- 
ally  inelastic  transitions  are  labelled  -*■  Jf.  Note 
that  higher  inelastic  channels  appear  whenever  the 
normal  component  of  the  incident  beam  energy  exceeds 
the  rotational  transition  energy. 


Figure  20  Rotational  transition  probabilities  PCJ  *  0  -►  J') 

J'  *  0,1, 2,3  as  a  function  of  E  for  HD,  calculated 

£ 

using  the  VEP  based  potential  with  anisotropy  para¬ 
meters  6^  =  +0.06,  Bg-  +0.09.  Potential  parameters 
for  VqqCz)  ate  given  in  Table  2.  Experimental  transi¬ 
tion  probabilities  are  indicated  by  crosses  (x)  . 
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Figure  21  Rotational  transition  probabilities  P(J  =  0  ■+  J') 

J'  =  0,1, 2, 3,  as  a  function  of  Ez  for  HD,  calculated 
using  the  exponential- 3  based  potential  with  aniso¬ 
tropy  parameters  8A  =  +0.06,  6R  =  +0.08.  Potential 
parameters  for  vj^Cz)  are  given  in  Table  2.  Experi¬ 
mental  transition  probabilities  are  indicated  by 
crosses  (x)  . 

Figure  22  Rotational  transition  probabilities  P(J  =  0  -+  J') 

J'  =  0,1, 2, 3  as  a  function  of  Ez  for  HD,  calculated 
using  the  exponential-3  based  potential  with  aniso¬ 
tropy  parameters  $A  *  +0.05,  BR  *  +0.20.  Potential 
parameters  for  Vqq(z)  are  given  in  Table  2.  Experi¬ 
mental  transition  probabilities  are  indicated  by 
crosses  (x) . 
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